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The large variations of T c across the cuprate families is one of the major unsolved puzzles in 
condensed matter physics, and is poorly understood. Although there appears to be a great deal 
of universality in the cuprates, several orders of magnitude changes in T c can be achieved through 
changes in the chemical composition and structure of the unit cell. In this paper we formulate 
a systematic examination of the variations in electron-phonon coupling to oxygen phonons in the 
cuprates, incorporating a number of effects arising from several aspects of chemical composition and 
doping across cuprate families. It is argued that the electron-phonon coupling is a very sensitive 
probe of the material-dependent variations of chemical structure, affecting the orbital character of 
the band crossing the Fermi level, the strength of local electric fields arising from structural-induced 
symmetry breaking, doping dependent changes in the underlying band structure, and ionicity of 
the crystal governing the ability of the material to screen c-axis perturbations. Using electrostatic 
Ewald calculations and known experimental structural data, we establish a connection between 
the material's maximal T c at optimal doping and the strength of coupling to c-axis modes. We 
demonstrate that materials with the largest coupling to the out-of-phase bond-buckling ("Bi 9 ") 
oxygen phonon branch also have the largest T c 's. In light of this observation we present model T c 
calculations using a two- well model where phonons work in conjunction with a dominant pairing 
interaction, presumably due to spin fluctuations, indicating how phonons can generate sizeable 
enhancements to T c despite the relatively small coupling strengths. Combined, these results can 
provide a natural framework for understanding the doping and material dependence of T c across 
the cuprates. 

PACS numbers: 74.72. Gh, 71.38.-k, 74.78.-w 
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I. INTRODUCTION 



Due to the extensive studies on the physical proper- 
ties of the cuprates, many constraints on the pairing 
mechanism of their high-temperature superconductivity 
(HTSC) have been accumulated. There is no doubt 
that the strong Coulomb interaction and the resultant 
strong electron correlations play crucial roles. This ef- 
fect is believed to be described by single-band Hub- 
bard or t-J models in 2D, and the magnetic mechanism 
for superconductivity has been proposed with the focus 
on the short range antifcrromagnctism or spin singlet 
formation. 1 These models have achieved great success in 
explaining many of the physical properties, such as the 
pseudogap, generalized magnetic susceptibility observed 
by neutron scattering, and the single particle Green's 
function found in angle-resolved photoemission spectra 
(ARPES). However, these models are not successful in 
explaining the variation of in superconducting transition 
temperature T c from material to material and other ma- 



terial dependent properties. For example, the famous T- 
linear resistivity within the plane is universally observed 
among various cuprates, while T c 's differ by two orders of 
magnitude. 2 ' 3 It is surprising that the resistivity, which 
is one of the most representative physical observables of 
the electronic states in solids, is irrelevant to T c . The 
2D Hubbard and t-J models contain only a few param- 
eters, such as hopping parameters t, t' , t" and interac- 
tions U and J. One possibility is that the range of the 
hopping and magnitude of t' and t" are key factors de- 
termining T c , which is determined by the structure and 
chemical composition perpendicular to the Cu02 plane. 4 
However, studies on the t-J model have found that fi- 
nite if suppresses superconducting correlations. 5 Recent 
investigations of the single-band Hubbard model using 
cluster dynamical mean field theory calculations 6 do not 
show increased tendencies towards pairing for larger if 
but variational studies do, 7 although the latter may be 
less controlled. 

From a structural point of view, the only known em- 
pirical rule of T c is that it increases at optimal doping 



as the number of CuC>2 layers n is increased for n < 3. 
Anderson noticed this n-dependence at a early stage, and 
proposed the interlayer mechanism of superconductivity 8 
The idea is that the single-particle interlayer hopping is 
suppressed by strong correlations within the layer, while 
the two-particle hopping is not. The onset of the latter 
below T c leads to the condensation energy of supercon- 
ductivity. Experimentally it is found that the c-axis elec- 
tron hopping is actually suppressed and there is no coher- 
ent band formation perpendicular to the plane. 3 There 
is no plasmon observed in the normal state, and below 
T c the Josephson plasmon appears in the low energy re- 
gion (~ 10 meV). 9 The idea of an interlayer mechanism 
has also been criticized in light of the c-axis oscillator 
strength, and has subsequently been abandoned. 10 

Increasing evidence for the importance of out-of-plane 
effects continues to accrue, which points to limitations 
of intrinsic planar models for the cuprates. For exam- 
ple, the correlation of out-of-phase oxygen dopant ions 
in Bi 2 Sr 2 CaCu208 + ,5 (Bi-2212) with features in the tun- 
neling density of states 11 has been interpreted in terms of 
a local increase of the superconducting pair potential. 12 
In addition, the rapid suppression of T c with out-of-plane 
cation dopants compared to in-plane dopants is surpris- 
ing given that the former do not appreciably affect in- 
plane resistivities. 13,14 T c has also been empirically cor- 
related with the Madelung energy difference between api- 
cal and planar oxygen atoms. 15 This has been recently 
supported by studies on Ba2Ca3Cu40sF2, a compound 
which has vastly different transition temperatures by ex- 
changing F with O at the apical site. 16 ARPES stud- 
ies have inferred a pairing gap which is a factor of 2 
larger on the bonding band in comparison to the anti- 
bonding band, 17 a trend consistent with other cuprates. 18 
Since the Fermi surface of the bonding band lies far away 
from either the antiferromagnetic reciprocal lattice zone 
boundary, or the van Hove points, linking the pairing 
mechanism with a purely electronic mechanism is not 
straightforward. 

As discussed in Refs. 4 and 15, the Madelung en- 
ergy difference and t' are directly linked. Pavarini et 
al. pointed out that the maximal T c in each family of 
cuprate materials scales with the next nearest hopping 
t'. 4 The energy of the Cu 4s orbital relative to the pd- 
a* band largely determines t' as the hybridization with 
the planar oxygen orbitals allows electrons to more effec- 
tively hop between 2p x>y orbitals. Strong apical 2p z -4s 
hybridization, determined by the Madelung energy differ- 
ence, raises the energy of the 4s orbital relative to the pd 
band, reducing the effective hopping t' . Thus, the further 
away the apical oxygen is located from the CuC>2 plane 
the larger i'. 4 Similar empirical relations between T c and 
structural details can also be found. For example, there 
is an optimal distance between the apical oxygen site and 
the mirror plane of the unit cell for which T c takes on its 
maximum value, as shown in Fig. 1. However, despite 
these observations, a direct connection between changes 
in the bandstructure and the pairing mechanism is lack- 




0.2 0.25 0.3 0.35 0.4 

Apical Distance [Units of c] 



FIG. 1: (Color online) The T c at optimal doping for many 
cuprate superconductors plotted as a function of the distance 
between the apical oxygen site and the mirror plane located 
at the center of the unit cell. The Hg- (red) circles, Tl- (blue) 
squares and Bi-families (green) triangles are shown as well as 
some LSCO and LBCO systems open (black). 



ing and to date, the connection between T c and these 
c-axis effects remains an empirical intrigue lacking a firm 
microscopic understanding. 

A phonon mechanism has also been pursued since the 
discovery of HTSC. 19 The original intuition by Bednorz 
and Miiller was that the Jahn- Teller (JT) effect with 
high-frequency oxygen phonon modes leads to HTSC. 20 
However, the degeneracy of the two e g orbitals is lifted 
considerably (~ 1 eV) and the JT effect is now be- 
lieved to be ineffective. Furthermore, there are several 
experiments suggesting a minor role played by phonons. 
These include the small isotope effect on T c at optimal 
doping, 21 the absence of the phonon effect on the tem- 
perature dependence of the resistivity, 2 and the absence 
of the phonon bottle-neck effect. 22 This conclusion was 
also supported by early density functional calculations 
which did not give appreciably high values of couplings 
and did not support a high T c in a BCS picture. 23 ' 24 

The high-T c cuprates, as strongly correlated systems, 
have attracted a great deal of theoretical interest. 25 ~ 41 
Part of this interest has been focused on examining the 
renormalization of electron-phonon (el-ph) interactions 
by strong electron-electron (el-el) correlations. 25 ~ 39 For 
example, studies of the t — J model incorporating an el- 
ph interaction indicate that polaron crossover may occur 
at a weaker el-ph coupling strength than in the case of 
a pure el-ph coupling model. 31 ~ 33 Quantum Monte Carlo 
treatments of the single-band Hubbard-Holstein model 
have shown that the renormalized el-ph vertex devel- 
ops a strong forward scattering peak with no substan- 
tial suppression of the el-ph vertex when the Hubbard 
U is large 34 with similar results obtained for the t-J 
model. 35 ' 36 However, using slave-boson approaches, the 
Hubbard model plus el-ph interaction was found to pos- 
sess no significant forward scattering peak together with 
an overall suppression of the el-ph vertex at low temper- 
ature. 37 This same study observed an enhancement in 



the el-ph vertex for small q scattering at high temper- 
ature, which was linked to phase separation. Cumulant 
expansion techniques 38 have also found an enhancement 
of the el-ph vertex for small q, again interpreted in terms 
of incipient phase separation when approaching a critical 
value of U. 

The single-band Hubbard-Holstein model has also been 
studied within dynamical mean field theory (DMFT). In 
one study 26 the effect of the el-ph interaction was to sta- 
bilize the insulating state in the vicinity of the density- 
driven Mott transition. In a paramagnetic DMFT 
study 28 the el-ph interaction was found to have little ef- 
fect on the low energy physics produced by the Hubbard 
interaction while modifying the spectral weight associ- 
ated with the upper and lower Hubbard bands at higher 
energy. From these results it was concluded that the pri- 
mary effect of the el-ph interaction was to reduce the 
effective value of U. Later work 29 considered antiferro- 
magnetic solutions in the presence of el-ph coupling and 
observed different behavior with strong polaronic effects. 
In this case, the critical coupling for polaron formation 
was observed to shift to larger values as the system was 
doped away from half-filling. The dichotomy between the 
paramagnetic and antiferromagnetic treatments 28,29 in- 
dicates the possible importance of the magnetic order in 
considering the el-ph interaction in correlated systems. 
Finally, a DMFT study invoking a Lang-Firsov trans- 
formation for the lattice degrees of freedom has found 
evidence for a competition between the Mott insulating 
(metallic) and bipolaronic insulating phases near (away 
from) half- filling. 30 

Using the dynamical cluster approximation, an ex- 
tension of DMFT, examinations of the el-ph interaction 
within small Hubbard clusters 39 find an overall suppres- 
sion of d-wave superconducting T c with increasing el-ph 
coupling. This occurs despite an increase in the apparent 
pairing correlations within the d x 2_ y 2 channel. The re- 
duction in T c was attributed to polaron formation, which 
reduces quasi-particle weight at the Fermi level and sup- 
presses T c through the loss of carrier mobility. The en- 
hancement of pairing correlations reported in Rcf. 39 in- 
dicates that the bare el-ph vertex has been renormalizcd 
in favour of (i-wave pairing, consistent with the obser- 
vations of Ref. 34. Furthermore, exact diagonalization 
(ED) studies on the t-J model, which include el-ph cou- 
pling to buckling and breathing vibrations, also show 
that the former enhance d-wave pairing while the latter 
suppress it. 25 While these results provide no definitive 
interpretation of the effect of el-ph coupling in strongly 
correlated systems, there is strong evidence that what- 
ever impact strong correlations may have, the el-ph in- 
teraction may still play a significant role that should not 
be overlooked in these systems. 

From an experimental front, the role of phonons in 
HTSC has become more prominent in recent years. 
While long-studied from Raman, infrared, and neu- 
tron measurements, recent ARPES and scanning tun- 
neling microscopy (STM) measurements on several 



cuprates has reinvigorated the exploration of the role 
of phonons on HTSC. 42 Experiments on Bi-2212, 43 ~ 46 
La 2 _ a; Sr a: Cu04 (LSCO), 47 Bi 2 Sr 2 Cu0 6+(5 (Bi-2201), 48 
Ba 2 Ca 3 Cu40 8 F2, 49 and Tl families Tl 2 Ba 2 CaCu20 8 (Tl- 
2212), TlBa 2 Ca 2 Cu 3 9 (Tl-1223), and Tl 2 Ba 2 Cu0 6 (Tl- 
2201) 50 have revealed kinks in the energy dispersion 
of these materials. These kinks have been interpreted 
as Hubbard renormalizations, 51 coupling to the neutron 
resonance and/or spin continuum, 52 and Engelsberg- 
Schrieffer renormalizations 53 due to coupling of electrons 
to a collection of optical phonons. These phonons in- 
clude the out-of-phase c-axis oxygen buckling modes and 
the in-plane Cu-0 bond-stretching modes. The disper- 
sion kink observed in the nodal region, (0, 0) - (n/a, it /a), 
and the peak-dip-hump structure observed in the anti- 
nodal region, (0, 7r/a)-(7r/a, 7r/a), clearly shows that the 
electrons are interacting with bosons of a well defined 
energy ~ 70 and ~ 36 meV, respectively. For the nodal 
region, it has been convincingly argued that this struc- 
ture is due to the oxygen bond-stretching phonon as the 
kink is observed independent of superconductivity and of 
the presence/absence of the spin resonance peak. As for 
the anti-nodal region, it has been claimed that the kink 
appears only below T c , and hence it is attributed to the 
spin resonance mode at 41 meV. 54 However, an exten- 
sive study conducted more recently has found the evi- 
dence for the kink structure in the normal state and over 
a wide range of momentum space. 44 Further, contrast- 
ing single and multilayer cuprate "kinks" , and materials 
known to have a neutron resonance, also indicates that 
the observed renormalizations are most likely due to op- 
tical phonons, although this is still controversial. 50 ' 55,56 

It is well known that the c-axis phonons show some 
of the most dramatic lineshapc changes with doping 
and temperature compared to any phonons observed via 
neutron 57 and Raman 58 scattering. For example, the 
apical phonon frequency shifts by as much as 20 cm -1 
with doping and temperature in a number of compounds: 
La 2 _ x Sr x Cu0 4 , HgBa 2 Ca„_iCu„0 4n+ 5 (n = 1 - 4) and 
Bi 2 Sr 2 Ca 2 Cu 3 Oio+,5. 57,58 Moreover, recent ARPES data 
on Bi-2201 have shown kinks in the energy range of the c- 
axis phonons which are weaker in overdoped compounds 
in comparison to optimal doped compounds. 48 This has 
also been interpreted in terms of increased screening of 
the el-ph interaction with increasing hole concentrations. 
In addition, the anomaly of the Raman Ai g -polarized 
mode due to the onset of superconductivity is observed 
in three- and four-layer compounds. 59 This has been suc- 
cessfully analyzed in terms of the internal electric field 
produced by the interlayer Josephson plasmon and its 
coupling to the phonon. This means that the system be- 
haves as an ionic crystal along the c-axis in the normal 
state, and suddenly turns into a superconductor. 

While the role of the neutron resonance and phonons 
remains controversial, it is of relevant interest whether 
these signatures in ARPES may be used as an angle- 
resolved analogy to the tunnelling ripples in conventional 
superconductors, 60 thus providing information on the 



pairing mechanism in the cuprates. In order to connect 
el-ph coupling to a possible pairing mechanism a system- 
atic study of coupling across families of cuprate materi- 
als is desirable. In Ref. 61 ARPES observed renormal- 
izations of the band were interpreted as due to the B\ g 
branch for anti-nodal electrons and the bond-stretching 
branch for nodal electrons. While the latter coupling 
is of a deformation type, the coupling constructed for 
the B\g branch involves a charge-transfer between pla- 
nar oxygens due to a modulation of the electrostatic 
or Madelung energies of the planar oxygen sites. A lo- 
cal crystal field, generated by a mirror plane symmetry 
breaking, allows for a coupling at first order in atomic 
displacements. 62 (We note here that the Ai g /Bi g nomen- 
clature only holds for Raman q — momentum transfers. 
However, throughout this work we denote the entire out- 
of-phase branch as u Bi g " and the entire in-phase branch 
as u Ai g " .) Since the cuprates are poor conductors along 
the c-axis, the electrostatic interaction can be thought 
to be largely unscreened. Calculations based on Ewald's 
method have been performed on YBa2Cu307 (YBCO) 63 
and large crystal fields have been obtained and the result- 
ing coupling matches well with the coupling determined 
from Fano lineshape analysis of Raman data. 62 

Recently, the issue of whether the el-ph coupling in 
the cuprates is strong enough to explain the observed 
band renormalizations has been revisited via density 
functional (LDA) calculations, 64 ~ 66 updating previous 
estimates. 23,24 While many efforts have been made to ex- 
tract bosonic coupling from ARPES renormalizations in 
the cuprates, there is no widely accepted way to uniquely 
determine the strength of the coupling at kink energies, 
and thus reliable comparisons of calculations with ex- 
periment must be viewed with some caution. 67 While 
LDA calculations have provided remarkably good agree- 
ment with phonon dispersions, there are a number of 
facets of LDA calculations which may only provide part 
of the story of el-ph coupling. As LDA calculations 
overestimate the itinerancy of the electrons, they de- 
scribe the cuprates as good metals even at half-filling. 
Moreover, the obtained interlayer transfer integral and 
hence the c-axis plasmon frequency is both coherent and 
much larger than the experimental observation. Both of 
these factors serve to overestimate the screening ability 
of the cuprates, especially in the undcrdoped region, and 
thus underestimates the strength of the el-ph interaction. 
This may be one of the reasons why LDA predicts smaller 
linewidths for the half-breathing oxygen bond-stretching 
modes and the apical oxygen modes, sometimes by more 
than one order of magnitude. 68 It is therefore not clear if 
these findings indicate that the el-ph coupling is small or 
that DFT-based approaches alone are inadequate for de- 
scribing the physics of the cuprates. As has been found in 
STM experiments, 11 a nano-scale inhomogeneous struc- 
ture on a length scale of 15 A exists universally in Bi-2212 
and YBCO. This length scale cannot be larger than the 
screening length and we can conclude that the screening 
length within the Cu02 plane is not shorter than 15 A, 



much longer than the Thomas-Fermi screening length of 
the typical metal (~ 0.5 A). Furthermore, charge transfer 
between the layers is almost prohibited. As a result of 
the transfer integral between the layer is proportional to 
[cos(k x a) — cos(k y a)] 2 , the opening of the pseudogap in 
the (n, 0) and (0, 7r) regions strongly suppresses the inter- 
layer hopping. Considering these discrepancies between 
LDA and experiments, one can imagine that the el-ph 
coupling is in reality much stronger than LDA predicts. 
Focusing on the buckling modes, theoretical considera- 
tion have been limited to the two-dimensional plane and 
the inter-layer Coulomb interaction has been neglected. 
This is usually justified in the metal since the screening 
length is much shorter than the interlayer distance. This 
is not the case in the cuprates. 

In this paper we provide a comprehensive and self- 
contained story on el-ph coupling to oxygen phonons as 
a function of doping across the cuprate families. We for- 
mulate a theory for el-ph coupling in the cuprates tak- 
ing into account the local environment around the Cu02 
planes, the poor screening of charge fluctuations out of 
the plane, doping-dependent band character variations, 
and structural differences across the cuprate supercon- 
ductors. 

The organization of this paper is as follows. In section 
II we present a general discussion of el-ph coupling to c- 
axis oxygen phonons. After summarizing prior work on 
the out-of-plane planar oxygen modes, we then provide 
a derivation of the coupling to modes involving c-axis 
apical oxygen motion. In section III we then discuss the 
anisotropy of the bare couplings and examine the total 
strength of the bare el-ph coupling and its contribution 
to the single-particle self-energy and <i-wave anomalous 
self-energy 

In section IV we develop the formalism for poor screen- 
ing and examine its implications for the anisotropy and 
overall magnitude of the rcnormalizcd el-ph vertices. Due 
to the poor c-axis conductivity we find that the c-axis 
phonons cannot be effectively screened for small in-plane 
momentum transfer q 2 £>- This effect becomes more pro- 
nounced as the crystal becomes more ionic and screening 
becomes increasingly inoperable in the underdoped side 
of the phase diagram. However, in the case of the B\ g 
modes the coupling is anomalously anti-screened produc- 
ing an enhancement of the coupling in the anti-nodal re- 
gion. In terms the projected d-wave couplings, the small 
q 2 £> behavior of the screened vertices produces an en- 
hancement in the total phonon contribution to pairing. 
Therefore, the combined effects of poor screening reduces 
the total el-ph coupling and enhances the d-wave pro- 
jected coupling with doping. This has important implica- 
tions for the doping dependence of the el-ph self-energies 
probed by ARPES as well as any contribution to pairing 
mediated by phonons. 

In section V we turn to materials trends and a sys- 
tematic examination of the Madelung potential and crys- 
tal field strengths across the Bi-, Tl- and Hg families of 
cuprates is presented. Here, using an ionic point charge 



model and the Ewald summation technique, we identify 
systematic trends in the strength of the crystal fields 
which mirror trends in the material's T c at optimal dop- 
ing. Through this observation we link the structure and 
chemical composition to the strength of the coupling to 
the c-axis modes and discuss how this can be used to 
understand the large variations in T c observed across 
the cuprates. We also present considerations for doping- 
induced changes to the value of the crystal field in Bi- 
2212. 

In light of these findings, section VI presents a simple 
two-channel model for pairing in the cuprates, which in- 
cludes a dominant, d-wave pairing, high-energy bosonic 
mode and a weaker phonon mode. Using this model, 
we demonstrate that phonons can provide a sizeable en- 
hancement to T c which is in excess of the T c that would 
be obtained from phonons alone. Furthermore, due to 
the dominant bosonic mode, the resulting value of the 
isotope exponent a is small (a < 0.15) despite the large 
enhancement of T c (~ 40 K). This calculation, in combi- 
nation with the materials and doping dependent trends 
identified in the previous sections, shows that a phonon 
assisted pairing model provides a natural framework for 
understanding trends observed for T c across the cuprates. 
Finally, in section VII we conclude by summarizing our 
findings and discuss open questions concerning el-ph cou- 
pling in strongly correlated systems. 

In addition to the treatment outline above, in the ap- 
pendix we explore how el-ph coupling to c-axis modes 
is modified by strong correlations in the half-filled par- 
ent insulators using exact diagonalization of small multi- 
band Hubbard clusters. Specifically we address how el-ph 
coupling modifies the properties of the Zhang-Rice sin- 
glet (ZRS). Here we find that static lattice displacements 
have a strong influence on the ZRS hoppings, energy and 
antiferromagnetic exchange energy J. These results have 
a direct impact on the use of down-folded models such 
as the t-J model as they indicate that the effects of the 
el-ph coupling cannot be simply cast as modulations of a 
single parameter such as t, J, or the energy of the ZRS. 



II. GENERAL ELECTRON-PHONON 
CONSIDERATIONS IN THE CUPRATES 

In this section we present a review of some generic con- 
siderations for electrons coupling to oxygen motions in- 
and out- of the CuC>2 plane. Since many of the deriva- 
tions have appeared before, we can be brief, with the 
main aim to generalize previous results to a five band 
model in order to include off-axis orbitals and apical oxy- 
gen phonon modes. 

We begin by considering an ideal CuC>2 plane isolated 
from its environment. Since hopping integrals arc modu- 
lated to second order in atomic displacements along the 
c-axis, the el-ph coupling due to this mechanism is weak. 
However, if the same plane is placed in an asymmet- 
ric electrostatic environment a local crystal field, which 



breaks mirror plane symmetry, provides a coupling lin- 
ear in displacement. The plane must then spontaneously 
buckle in a pattern where the oxygen (copper) atoms are 
displaced away from (towards) the mirror plane. 

In a three band model, 62 the local field coupling was 
used to construct a charge-transfer el-ph vertex for cou- 
pling to the Raman active out-of-phase and in-phase c- 
axis oxygen vibrations. The in-phase phonon modulates 
charge transfer between planar oxygen and copper or- 
bitals while the out-of-phase phonon modulates charge 
transfer between only the planar oxygen orbitals. A 
single-band cl-ph coupling was obtained: 
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| ff (k,q)| 2 c3 c _ q , CT c k , CT [6t 1 + 6_ q ], (1) 



where c k (ck,cr) creates (annihilates) an electron in the 
partially filled antibonding band with momentum k, en- 
ergy £k) and spin a, and 6 q (6 q ) creates (annihilates) a 
phonon of energy 17 q and wavevector q. 

Considering a modulation of the electrostatic coupling 
of the charge density at the oxygen sites coupled to the 
local on-site potential & e xt, the Hamiltonian is of the 
form 62 
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AArf[u(an)], 
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where u(an) is the oxygen displacement vector in the 
unit cell at lattice site n, e is the electron charge and p^ n a 
(Pn,<r) is the creates (annihilates) an electron at site n, 
which can include both planar and apical oxygen orbitals. 
This coupling mechanism differs from the deformation 
coupling considered in Ref. 24. Expanding for small 
displacements H' site = H sUe +H e i- ph +0(u 2 ) where H site 
includes the Madelung contribution to the site energies 
and the term linear in u generates the el-ph interaction 
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(3) 



E n is the local crystal field at the oxygen site provided 
this field is finite, which occurs at locations of broken 
mirror symmetry in the unit cell. 

In order to derive the form of the coupling g(k, q) Eq. 
(3) must be rewritten in the form of Eq. (1). To do so, the 
el-ph Hamiltonian is Fourier transformed to momentum 
space and the oxygen operators are replaced by band 
representation operators p^s.a- — 05(k)ck, CT - Here, 0,5 (k) 
is the oxygen (<5 = x, y for planar oxygen and 5 = a for 
apical oxygen) eigenfunction for the pd-cr* band, which is 
obtained from a tight-binding model for the CL1O2 plane. 



A. Multi-band Models 

Prior work focusing on the A\ g , B\ g and breathing 
branches made use of a three-band model. 15 ' 61,62 In or- 
der to extend these works to include the apical oxygen 




FIG. 2: (Color online) The five-band model used to derive 
the form of the el-ph couplings g(k, q). 



modes this model must be extended to a five-band model 
as shown in Fig. 2. The basis set of this model con- 
tains a 4s (s n)CT , sj, a ) and 3d x 2_ y 2 (d a ^ a ,d^ n a ) orbital 
on each copper site n, two planar oxygen 2p x ^ y orbitals 
(Pnj^^g^) with 5 = x,y, and one apical oxygen 2p z 
orbital (an^aj, CT ). Here, we neglect the in-plane O 2p or- 
bitals oriented perpendicular to the Cu-0 bonds. These 
orbitals form weaker pd-n bonds with the lower energy 
Cu tig orbitals and do not contribute heavily to the char- 
acter of the band crossing the fermi level. 24 Site energies 
are denoted by e s ^, P ,z, respectively. Defining canonical 
fermions 69 a, (3 from combinations of the planar oxygen 
orbitals via a Wannier transformation 



n ,.S x ,y{k-)Pis.,x,(yTSy } x(>i)Pis.,y,iT , ,n 

ak,<T,Pk,<T = ±% , (4j 

Mk 

where s XiV = sin(k x . y a/2) and /x 2 . = s 2 (k) + s 2 (k), the 
Hamiltonian can be written as H = ^ k iJ k a : 



-Hk,<7 = Hsite - 2tpj ) tl c [/3k, -/3k,a- - Oi^OC^g-] 

+ 2 



t P dHk,<rdl k ff a k , ff + tppXkalefa., 



t' 



p s J 



i P s K kSk ](T ak, CT - *p S AkSk i0 ./3k, CT + -^- s k, CT a k,o- 



tpp^at a a K<? + t'pp^aL aPk,a + h.C 



with H S it e containing the site energies and h.c. denoting 
the hermitian conjugate. Finally, the basis functions are 
defined as: 



v*. 



,4(k)4(k) 



M k 



K k 



4(k)-^(k) 



/'-k 



A k = 2 



s x (k)s y (k) 



/'k 



Xk = A k «k- 



(6) 



In the limit where the apical and copper 4s orbitals are 
removed (t ps — t' = 0, etc.) there expressions recover 
prior work carried out using a three-band model. 15 ' 61 ' 6 



B. Planar Oxygen c-Axis Modes 

For planar oxygen vibrations the el-ph coupling is given 
by (p = k - q) 



9Bi g ,Aig(k,q) = eE 



^ 2M N((i)n BlgAlg 
4(k)&.(p)e -, ' 9 " a/a (l + e- 1 ' 9 *«)=F 

<^(k)^(p) e -^W 2 (l + e-^ a ) 



(7) 



where the minus (plus) sign is for coupling to the Bi g 
(Ai g ) branch. Here Mo denotes the oxygen mass, E z 
is the c-axis component of the local field at the planar 
oxygen site, A^ 2 (q) = 4 [cos 2 (q x a/2) + cos 2 (q y a/2)] is the 
phonon eigenvector normalization, and ^Aig,Big denote 
the assumed dispersionless frequencies of the B\ g and 
A\g branches, respectively. The motion of the heavier 
Cu atoms has been neglected in this treatment. Eq. (7) 
is a generic expression for the A\ g and B lg vertices, in- 
dependent of the underlying tightbinding model used to 
determine the band eigenfunctions <^ x ,j/,Cu(k). In the 
three-band model, the band eigenfunctions are defined 



,61 



Wk) = Tj^T[e(k)t x , y (k)-t'(k)t v , x (k)] 
^cu(k) = -^[e 2 (k)-i' 2 (k)], 



(8) 



where t x>y (k) = t pd s x , y (k), t'(k) = -4i w ,sin a .(k)sin 1/ (k), 
the normalization is 



A 2 (k) = [e 2 (k)~t' 2 (k)] 2 + [e(k)t x (k) 
+ [c(k)*„(k)-t / (k)t a .(k)] 2 1 



t'(k)i,(k)] 5 



(9) 



and e(k) is the bare dispersion, given in Ref. 62. From 
Eq. (7) it can be seen that the symmetry of the phonon 
is implanted into the el-ph coupling to provide substan- 
tial momentum anisotropy. For the B\ g phonon branch, 
g(k, q) changes sign for k x ,q x — > k y , q y , while for the A\ g 
phonon branch it does not. As a result the coupling in- 
(y) volves all fermionic states for the A ig branch, while for 
the Big branch the antinodal states along the Brillouin 
zone (BZ) axes are weighted heavily and nodal states 
along the zone diagonal are projected away. 



C. Apical Oxygen c-Axis Modes 

Since apical phonons show some of the strongest renor- 
malizations in La- and Hg-cuprates, 57,58 and since the 
apical oxygen atoms do not lie in a mirror plane sym- 
metry even in single layer cuprates, they are included as 
an extension of our previous work. Early on the apical 
phonon was thought to be quite anharmonic and related 



to the Jahn- Teller mechanism in YBa2Cu307, 70 although 
many confusing results were found. 71 More recently the 
coupling is thought to be electrostatic 72,73 in nature and 
in some treatments weakly momentum dependent. 72 It 
should be emphasized that formally a Holstein (momen- 
tum independent) coupling in any model with long-range 
Coulomb interactions will be screened out by backflow 
due to charge conservation and therefore the coupling is 
expected to be very small. This will be discussed in Sec- 
tion IV. Here instead we place focus on a strongly mo- 
mentum dependent coupling arising from charge trans- 
fer mechanisms between apical and planar oxygen or- 
bitals, similar to the electronic pathways involved in c- 
axis tunnelling. 74 

An apical orbital displacement modulates the 
Madelung energy as in Eq. (3) and the resulting el-ph 
coupling is of the form of Eq. (1), with 



ffa pe:c (k,q) = g p </>T(k)</> a (k-q)e*(q) 



(10) 



Here, e„(q) is the c-axis component of the eigenvector 
for the apical branch coming from the atomic displace- 
ment in Eq. (3), gl pex = eE% v / h/2M fl a , and o (k) 
is the apical eigenfunction, obtained from diagonalizing 
Eq. (5). For our purpose however, we are primarily in- 
terested in the leading order momentum dependence of 
the coupling g(k,q). Starting from Eq. (5), the Lowdin 



down-folding procedure 75 is applied to determine the api- 
cal character of the resulting partially filled band crossing 
the Fermi level. The resulting form for the apical eigen- 
function (p a (k) is then 



a (k) = 2t pz 



K k 



Ck 



with «k defined in Eq. (6). Finally, for simplicity we 
neglect the momentum dependence of the apical phonon's 
eigenvector and set e*(q) = 1. More complicated models 
in which the apical phonon involves the motion of the 
in-plane oxygens can be treated accordingly. 

The local crystal field E® at the apical oxygen site 
modulates a charge transfer between the apical oxygen 
and the planar orbitals. This mechanism of charge trans- 
fer is analogous to the charge transfer mechanism yielding 
bi-layer splitting. 74 From Eqs. (5) and (10), the resulting 
momentum dependence of the coupling via this transfer, 
g(k, q) ~ [cos(fc x a) — cos(fc y a)][cos(p 2: a) — cos^a)] with 
p = k — q, is strongest for anti-nodal electrons, and has a 
form factor similar to c-axis hopping £j_(k). 74 Although 
this coupling has an anisotropy similar to that of the cou- 
pling to the B\g branch, it does not contribute to <i-wave 
pairing due to its phase and momentum dependence at 
large q. 



D. In-plane Bond-Stretching Modes 



For completeness, we also consider the coupling to the planar Cu-0 bond stretching modes, the so-called breathing 
modes, within the framework of the three-band model, as derived in Refs. 32, 36 and 61. The bond-stretching modes 
couple to electrons via both a direct modulation of the hopping integral t p( i as well as electrostatic changes in the 
Madelung energies as the orbitals are displaced. 36 As done in Ref. 32 we consider only the overlap modulation. The 
derivation is briefly sketched here. 

To obtain the form of the el-ph coupling the overlap integral t p d is taken to be site dependent t pd . It is then assumed 

that the Cu and O atomic displacements, u^ u and u° al about their equilibrium positions, R n and R n + la/2, where 

S are basis vectors for the CuC>2 plane, are small. The overlap integral is expanded and only the first order term is 
retained 



pd ^pd 



l r,d T 



E 

=±x,±y 



Vt pd\ r =R, 



■ (u Cu 



ha<5/2 



0(u 2 



The modulation of the hopping integrals provides the el-ph coupling Hamiltonian 



H el-ph — X^i P^t pt 



pd | 



Cu 



1° - 

n+a8/2 



[4 



aPr, 



h.c. 



(11) 



(12) 



with P X]y = ±1 = —P-x-y denoting the phase of the Cu-0 overlap. Following Ref. 61, we neglect the Cu vibration 
and set dt^ d /dx\n n = —Qs9dp, where Q± x = Q± y = ±1 and g^ p is a scalar function that depends on the equilibrium 
Cu-0 distance. The el-ph coupling Hamiltonian can then be simplified to 



H el-ph = 9dp Yl P sQsuf( n ) [ d t,aPn,a,S + h.c] 



(13) 



as obtained in Ref. 61. Here v? s denotes the displacement of the oxygen atom S along the Cu-0 bond. By introducing 
the Fourier transform the el-ph coupling may be re- written as: 



h: 



br 
el—ph 



br 



^2 ^2 p s cos(k a a/2] 

k.q a,8— x,y 



tf(q)[4 



k.<jPk-q,(T,(5 



■A- 



q,(j,<5 



rfk,a](6l + &-q)- 



(14) 



Here, e^(q) = sin(qsa/2) / \i sin (q x a/2) + sin (q y a/2) is the component of the phonon eigenvector for oxygen (5 parallel 

to the Cu-0 bond. Finally, the electronic eigenfunctions of the planar O (f> x . y (k) and Cu 3d 2 — y 2 orbitals </>cu(k) for 
the pd-cr* band, with operators c, c\ are introduced. The resulting Hamiltonian reduces to the form of Eq. (3) with 



5br(k,q) = 3o r Yl ^a^W [ C0S iP a a/2)(f> Cu (p)(j) a (k) - cos(fc Q a/2)0 c „(k)</> Q (p) 



(15) 



a—x.y 



9 h J 



= gd P y // h/2MoQq- We note that gen- 



Here p = k— q, 

erally t pd (d Cu -o) ~ d P Cu -o and tnus 5 P d = Pt p d/d Cu -o- 
Since typically /3 = 3.5, 32 an estimate for the strength of 
the coupling is obtained from g pd ~ 2 cV/A for t p d = 1.1 
eV. For tt„ = 70 meV, this gives gt, r — 86 meV. 



III. K,Q-MOMENTUM DEPENDENCE OF THE 
BARE VERTICES 



A. Momentum dependence throughout the 
Brillouin Zone 



In the previous section it was shown how the explicit 
form for el-ph coupling to oxygen modes is determined 
by the nature of the charge transfer modulated by the 
lattice displacement, the local environment surrounding 
the Cu02 plane, as well as the orbital content of a single 
downfoldcd band crossing the Fermi level. The relevant 
parameters - magnitude of orbital hybridization, local 
crystal field, the charge-transfer energy, the shape of the 
Fermi surface, and the density of states at the Fermi 
level - all contribute in setting the overall magnitude of 
the coupling as well as the full fermionic k and bosonic 
q momentum dependence of the coupling g(k, q). The 
band character enters through the band eigenvectors 4> 
which further depend on the complexity of the unit cell. 
At the BZ center the wavefunctions are atomic and the 
character of the band is unique. Large momentum vari- 
ations of the band character then occur for increasing 
momentum, and a very strong momentum dependence 
of the overall el-ph coupling can occur. This strong mo- 
mentum dependence has indeed been observed in recent 
LDA treatments 64 . 

We remark that we are first interested in the magni- 
tude and anisotropy of the bare couplings in the absence 
of charge screening in order to determine possible discrep- 
ancies with LDA treatments, which treat correlations on 
the mean field level and give 3D metallic screening. In 
order to estimate general tendencies, momentum depen- 
dencies as well as magnitudes, in this section we explore 
some simplifications. 



We begin by assuming that 2t pc i is much greater than 
any relevant energy scale in the system, keeping in mind 
that the charge transfer energy A = e p — e d =~ 0.8 eV is 
much reduced from its bare value ~ 3.5 eV when treating 
correlations in mean field approaches like LDA. 24 In this 
limit the <f> functions can be represented as 



</>x,3/(k) = ±i 



and 



sin(k Xiy a/2) 



sin 2 (k x a/2) + sin 2 (k y a/2) 



■ 0{A/2t pd ) 



OCu 



A/2t pd 



Jsm 2 (k x a/2) + sm 2 (k y a/2) 



To the same order, the denominators in these expres- 
sions are constant over constant energy contours. There- 
fore, since we will restrict ourselves largely to the Fermi 
surface, we represent the band functions (j) x _ y (k) — 
A sin(fc Q a/2), b (k) = A Cu and o (k) = A a d k , 
with dk — [cos(k x a/2) — cos(k y a/2)]/2 and the coeffi- 

(<t>l y (k))/(sm 2 (k x ^ y a/2)), 



cients determined by A 2 Q 



4 2 



b 2 (k)) and A\ = (^(k))/(d 2 ), respectively. 



Here (...) denotes a Fermi surface average: (A) = 
Sk^ k ^( £k )/Sk ^( e k)- In this way, the overall coupling 
anisotropy can be simplified without loss of generality 

As a consequence, the fermionic momentum depen- 
dence of the coupling to the breathing modes disappears: 

,9hr(k, q) = g h Q A Cu A ^ Pae a (q)sm(q a a/2). 



Substituting the phonon eigenvectors, the coupling to the 
breathing modes becomes 

3br(q) = g^AcuAo \Jsin 2 (q x a/2) + sm 2 (q y a/2). (16) 

This form has also been obtained in a t-J approach 32 
however, in this case, the oxygen and copper charac- 
ter have been explicitly retained through Ao and Ac u , 
respectively. 76 



The el-ph vertex for the A\ g and B\ g modes can be likewise simplified 



5Ai g ,sig(k,q) = eE, 



2h 



A 2 p-Hlx+ly) -/ 2 



M N(q)Sl BlgiAlg 

x [sin(fcj;a/2) sm(p x a/2) cos(q y a/2) ± sm(k y a/2) sm(p y a/2) cos(q x a/2)] 



(17) 



These expressions recover the Raman form factors in the limit q — > for each mode and they obey the symmetry 
conditions for momentum reflections about 45° as discussed previously. The fcrmionic momentum dependence cannot 
be neglected in either of these expressions, where in particular, a strong fcrmionic momentum dependence of the 
coupling to the B\ g modes occur, preferentially weighting anti-nodal states with small momentum transfers. 

Lastly, the momentum structure of the apical coupling simplifies considerably in the same manner: 



g ap ex(K q) = g" pex A 2 a [cos(k x a) - cos(fc y a)][cos(p £C a) - cos(p y a)}/4. 



(18) 



Once again, a substantial fcrmionic momentum dependence emerges from the c-axis charge transfer pathways and the 
apical character of the band. 



Before proceeding further a few comments are in or- 
der. In using Eqs.(16)-(18) we have simplified the cl-ph 
couplings while explicitly retaining the role of the band 
character in determining the overall strength of the cou- 
plings. Since the eigenfunctions enter to the fourth power 
for |g| 2 , the total coupling strengths determined in this 
approach may change considerably when adjusting multi- 
band parameters. However, this method has the advan- 
tage that the materials dependence of the coupling, pa- 
rameterized by Ao, A a and Ac u , can be calculated using 
a variety of methods such as exact diagonalization, quan- 
tum Monte Carlo or LDA. We also emphasize that the 
total coupling strengths (calculated in the next section) 
that we obtain using this formalism are similar to those 
obtained from LDA treatments, even though the latter 
includes the effects of screening. 23 > 64 > 65 

This approach also allows for the use of a renormalized 
bandstructure, while retaining the explicit band char- 
acter of the original five-band model. This is impor- 
tant since the overall strength of the el-ph couplings 
scales with the density of states at the Fermi level Np 
(see Eq. (19)); narrow bandwidth systems will exhibit 
larger coupling in comparison to large bandwidth sys- 
tems with the same vertex g(k, q). With an appropriate 
choice in parameters, the five-band model given in sec- 
tion II reasonably reproduces the bandwidth (and Np) 
determined by LDA calculations. 77 However, as has been 
noted, 68 LDA over predicts the total bandwidth (and 
consequently Np is under predicted) in comparison with 
experiment. Therefore, we expect that the total cou- 
plings will be underestimated using the five-band model 
with parameters chosen to match LDA. A simple rescal- 
ing of the five-band model bandstructure in conjunction 
with the full form of the (^-functions is insufficient to cor- 
rect this since this procedure would produce incorrect 
values for the (^-functions and therefore generate errors 
in the character of the band as a function of k. However, 
the use of the couplings defined by Eqs. (16)-(18) allows 
us to resolve this issue. Here, the correct band character 
is captured by calculating Aq , Aq u and A a using the five- 



band model but the Fermi surface and bandstructure are 
obtained from elsewhere in order to better match experi- 
ment. In this work, we adopt a five-parameter tightbind- 
ing model for Bi-2212 derived from fits to ARPES data. 78 
This approach allows us to capture the increased value 
of Np while simultaneously retaining estimates for the 
correct band character. We also note that the specific 
shape of the Fermi surface is not crucial to the overall 
anisotropy of the couplings. 



To visualize the momentum dependence of the cou- 
pling in more detail, we plot in Fig. 3 g(k-F, q), given by 
Eqs. (16) - (18), as a function of transferred momenta 
q along two directions as indicated. The dependency 
on transferred momenta arises from the nature of the 
charge-transfer coupling of the different modes. The c- 
axis modes (Figs 3al-2,4 and 3bl-2,4), being electrostatic 
in nature, translate into stronger coupling for small mo- 
mentum transfers, while the deformation-type coupling 
of the breathing branches gives stronger coupling at large 
q, and vanishes in the limit q — > (Figs. 3a3 and 3b3). 



Apart from the breathing modes, an appreciable 
fcrmionic wavevector dependence of the couplings to the 
c-axis modes is found as a consequence of the charac- 
ter of the underlying atomic vibrations. For the case of 
the A\ g and apical modes, for momentum transfers along 
the zone diagonal (Figs. 3bl and 3b4, respectively), the 
fcrmionic dependence is symmetric with respect to re- 
flections about 7r/4 while the B\ g coupling (Fig. 3b2) 
changes sign. Momentum transfers along the zone face 
(Figs. 3al, 3a2, 3a4) do not obey any set selection rule, 
although the symmetric- or anti-symmctric-like charac- 
ter of the coupling is evident. Finally, the strong momen- 
tum dependence of the charge transfer along the c-axis 
dictates that the apical coupling vanishes for any fermion 
momentum along the zone diagonal (Fig. 3b4). 
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FIG. 3: (Color online) Plots of the el-ph coupling constant 
g(ki?,q) for fermionic momentum on the Fermi surface as a 
function of transferred momentum (q x ,q v — 0) (al-a4) and 
q x = q y (bl-b4), respectively. al,bl (a2,b2) plot coupling to 
the Ai g (Big) branches, respectively, a3, b3 plot coupling to 
the breathing branch and a4,b4, plots coupling to the apical 
branch. The colors denote angles from the corner of the BZ 
(shown in the inset of Fig. 4) given by solid black - 0° , solid 
red - 15°, solid green - 30°, solid blue - 45° (nodal), dashed 
green - 60°, dashed red - 75° and dashed black - 90° (anti- 
nodal). 



B. Momentum dependence on the Fermi surface 

The strong dependence of g(k, q) on both k and q leads 
to anisotropic coupling between electrons and phonons. 
As the most relevant scattering processes involve those 
states near the Fermi level, the explicit momentum de- 
pendence as seen in ARPES is most clearly envisioned 
by calculating fc-dependent self-energies, in terms of el- 
ph coupling A. We take the modes to be dispersionless 



a 



l^v 



Vt Ul which couple to electrons via the Fock piece 



of the electron self-energy A^(k) in lowest order 79 : 



A„(k) = 



MX 



5>(k,q)| 2 ,5(£ p ) 



^<l*(k,q)l 9 >P„ 



(19) 



with £ p = e(p) — fi, fi is the chemical potential, N is the 
number of momentum points and Np = -4 X)k^(^k) is 
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FIG. 4: (Color online) Plots of A y (k_F) for momentum points 
along the Fermi surface. Inset: the definition of the angle 9 
used in Figs. 3 and 4. The parameters used are defined in 
the text. 



the density of states at the Fermi level. The delta func- 
tion restricts the sum to initial and final fcrmion states k, 
p, that lie on the Fermi surface, with scattering between 
them governed by the transferred phonon momentum q. 
The resulting A„(ki?) for the four modes are shown in 
Fig. 4. The Fermi surface contour and bandstructure 
are again determined from a five-parameter tightbinding 
bandstructure. 78 The el-ph coupling vertices are evalu- 
ated using Eqs. (16) - (18), along with the conventional 
parameter set (in eV): t p d — 1.1, t pp = 0.5, t pz = 0.29, 



Sp 



2, t s 



1.5, <i d ~ P = 0.8, e d - z = 1.00, e d _ s = 



_7_i5,80 p or these parameters we obtain Aq = 0.446, 
A 2 Cu = 0.592 and A\ = 4.52 x 10~ 2 . Additionally, esti- 
mates of the local field strength at the planar and api- 
cal oxygen sites for Bi-2212 are used: Ef ane = 3.56, 
jjjapex _ ig_33 eV/A, respectively. A systematic deriva- 
tions of these E field values across the cuprate families is 
given in Section V. 

A strongly varying A(k) is obtained for the c-axis 
modes, largely weighting anti-nodal fermion states, aris- 
ing from both the fermion dependence of the bare cou- 
plings as well as the small momentum transfers connect- 
ing anti-nodal points on the Fermi surface, as noted in 
prior treatments 61 as well as more recent LDA- projected 
studies. 64 The anisotropy is particularly strong for the 
apical coupling because of the strong fermionic momen- 
tum dependence of <?(k, q). On the other hand, coupling 
to the breathing modes has a much weaker anisotropy 
along the Fermi surface, due to large momentum trans- 
fers (7r, 7r) connecting anti- nodal portions of the Fermi 
surface, as well as (7r,0) and (0,7r) transfers connecting 
nodal points. Inclusion of a variation of the copper char- 
acter across the Fermi surface yields a more anisotropic 
coupling. 61 However, the character of the anisotropy is 
not as strong as that observed in Ref. 64, where the 
breathing branches couple more strongly to anti-nodal 



11 



Branch 


A, 


A^> 


Mg 


7.74 x 1(T 2 


4.42 x 10" 2 


Big 


0.17 


0.12 


Breathing 


8.27 x 10~ 2 


-4.75 x 10~ 2 


Apical 


l.n x i(r 2 





Total 


0.341 


0.1149 



TABLE I: Tabulated values of Az,,/, for the four phonon 
branches considered in this work. 



states. 

The full coupling A z , renormalizing the single particle 
self-energy, is given by a sum over all modes averaged over 
the Fermi surface, A z = X^A^k)), an< ^ can be visual- 
ized most easily as an average over the curves shown in 
Fig. 4. The renormalization A^ for a d x 2_ y 2--wave super- 
conductor gives the el-ph contribution to the anomalous 
self-energy and is given by a d-wave projected average 






£k,p<krfp|g,,(k,q)| 2 S(6 c )*(gp) 



(20) 



where q = p — k and d-wave basis function d^ — 
[cos(k x a/2) — cos(k y a/2)]/2, as before. A positive (nega- 
tive) A<£ denotes an attractive (repulsive) pair interaction 
in the d x 2_ y 2 channel. The A values obtained for the four 
phonon branches considered in this work are given in ta- 
ble I. As noted before, the A lg and B ig modes enhance 
d-w&ve pairing, the breathing modes suppress it and the 
apex modes gives no contribution. We note that these 
values for the X z are similar to those obtained from re- 
cent LDA results for YBa2Cu307 64 ' 65 but arc smaller 
than those calculated for doped CaCu02. 23 ' 24 We again 
remark that important consequences of screening are not 
considered here, which can affect both the magnitudes 
and anisotropies of the couplings as will be discussed in 
Sec. IV. 



C. Kinematic constraints 

It is important to note that in ARPES there is a strong 
kinematic constraint that also governs which fermion 
states can couple to each mode. Due to energy conserva- 
tion of scattering from dispersionless phonons, the nor- 
mal state ARPES spectra may show mixing of electron 
and phonon states if the band energy along a particular 
cut crosses the phonon mode energy. This gives rise to 
a kink in the dispersion, with trailing intensity observed 
asymptotically to the phonon mode energy. Otherwise, 
if the band along a particular cut lies above the mode 
energy, the spectra may show level repulsion in the form 
of a flattened band bottom along that cut. 

In the cuprates, the shallow band near the anti-nodal 
parts of the BZ limits the couplings of those states to 
only modes with energies less than £(0,7r). Thus, quite 
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FIG. 5: (Color online) A contour plot of the band structure 
of optimal doped Bi-2212 obtained from the five-parameter 
tight binding model of Ref. 78. The outermost contour (solid 
black) corresponds to the FS. The remaining contours corre- 
spond to £k — A* = Qv f° r the four phonon branches considered 
in this work. Reading from the FS to the T-point, the contours 
correspond to the B\ g , A\ g , breathing and apical branches, 
respectively. 



generally, we would expect coupling to the modes to dis- 
appear from the spectra when £7„ < £(k) along the entire 
cut. The contours £(k) = 0„, which define this kinematic 
constraint, are plotted in Fig. 5 at optimal doping. The 
Fermi level is also shown for reference and corresponds 
to the outermost contour (solid black). In the case of 
the apical mode, although the coupling may be largest 
for \s-an, a kink effect is prevented kincmatically. More- 
over, the breathing modes will be observable only for 
near-nodal cuts due to the form of the coupling as well 
as the energies of the modes which are large when com- 
pared to the A\g /Big modes and the saddle point energy 
£(0,7r). In fact, for optimally doped Bi-2212, the Bi g 
modes lie right at £(0, 7r), and can be seen most clearly 
in the anti-nodal directions, although its coupling should 
be felt across the entire Fermi surface as shown in Fig. 
4. 

In the superconducting state, the energy contours shift 
by the maximum value of the energy gap and are set by 
fl u + Ao, and the kinematic range shrinks for coupling to 
the higher energy modes. 61 ' 81 



IV. POOR SCREENING ALONG THE c-AXIS 

A. Overall Approach 

Screening - responsible for reducing the net el-ph cou- 
pling in conventional good metals - must be included for 
all c-axis phonons. In the usual treatment of screening, 82 
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the screened el-ph coupling g is 

V(q)U gtl (q,n) 



g(k, q, £2) = ,g(k, q) 



1-V(gf)ni,i(«,fi) 



(21) 



where V(q) = 4ire 2 /q 2 is the 3D Coulomb interaction and 
Yla t b{q, £2) is the frequency-dependent polarizability cal- 
culated with vertices a, 6, respectively. Since the plasmon 
frequency £2 p ; is usually much larger than the phonon fre- 
quencies in 3D metals, in the limit q — > the effective 
coupling is g(k,q — > 0) — Sg with <5<? denoting the av- 
erage value of g(k,q) over the Fermi surface. Thus the 
bare el-ph coupling is normally well screened. 

Turning now to the cuprates, from the point of view 
of Coulomb interactions, the materials are 3D with q 2 — 
q 2 D+q 2 and q 2D = (q Xl q y ) . However, due to the largely 
incoherent c-axis transport observed across the phase di- 
agram (apart from the overdoped side), the polarizabil- 
ity is largely determined by the planar conduction elec- 
trons and II ai f,(<7, £2) = n^(q 2 £>,£2)- As a result, impor- 
tant changes to the effects of screening occur for small 

q 2 c < q- 

A qualitative feel for the effects of poor screening can 
be obtained by considering a coupling which is inde- 
pendent of Fermion momentum k and in the limit of 
small in-plane momentum transfers, where VFq2D <C 
Q p h- This limit is most relevant for Raman active 
phonons and n 11 (g, £2 p / l ) = nq 2 D /mfl 2 hl with n/m = 

VKdiEk/faO^ and n ff ,i(g,£2 P fc) = -%k , with 



q2 

n/m g (q) = 2/V ce nEk/( e k)5(q)^ where V ce u is the 
unit cell volume and / is the Fermi distribution. The 
screened el-ph interaction is then 



s(k,q,£2 p /0 =g(q) 



m 



K(<i) 



m g(<l) n lM) 



<V 



(22) 



with £2pj(<7) = £2p ( • {qiD/q) 2 - Thus c-axis Raman-active 
phonons in the cuprates, with couplings strongest for 
small momentum transfers q and Q p i(q) < £2 p /j, should 
survive the effects of screening. 



B. Considerations for the polarizability and \ v 

To evaluate the polarizabilities in the cuprates, we use 
the standard Lindhard expression 



n a ,b(g,£2) 



v, 



cell 



5>(k,q)&(k,q) 



/(Ck) - /(£k+q) 



£k — e k+q 



+ £2 + ic5' 



(23) 



with all wavevectors being two-dimensional. This form of 
the polarizability has been considered in previous works 
focusing on the role of screening and reduced dimension- 
ality in the doping-dependent softening of the in-plane 
bond-stretching modes. 83,84 For the Coulomb interaction 
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FIG. 6: (Color online) (a) Plots of the real (black) and imag- 
inary (— lx, red) parts of the charge susceptibility ITi,i and 
(b)-(f) the renormalized el-ph coupling g(k, q) in arb. units 
for k = (0, 7r/a) and q = (q x ,0) for different values of q z 
as indicated by the legend. For (b)-(f) the black dashed 
line corresponds to the bare coupling g(k, q). A momentum- 
independent bare coupling is considered in (b) , followed by the 
Ai g , Big, apical and breathing modes in (c)-(f), respectively. 
In the case of the breathing modes a small q y component has 
been added to visualize the effects of screening at small q 2D . 
The q x axis in all plots is on a logarithmic scale to highlight 
the small q 2r) region. 



V(q), more appropriately for the cuprates, we consider a 
layered system with charge density essentially confined to 
2D layers, Coulombically coupled across the planes: 85 ' 86 



V(q a b,q z ) 



V 



1 



q ab ateLnh(q ab /q ) 1 + F 2 



(24) 



where V = 27re 2 /a v /e a6 e c , q Q = (4/c)- v /e c /e a6 with e abyC 
the in-plane, out-of-plane dielectric constants, respec- 
tively. The c-axis lattice constant sets the in-plane mo- 
mentum scale beyond which the planes become effec- 
tively uncoupled, and leads to an interpolation between 
the 2D and 3D form for the Coulomb interaction where 
F z = sm(q z c/A)/ smh(q a b/qo) provides the q z dispersion 
of the interaction. 

The real and imaginary parts of the pure charge po- 
larizability Il^i are plotted in Fig. 6a. Here, S = 5 meV 
has been included to smooth singularities, e(k) is taken 
for optimal doping in Bi-2212, 78 c = 30.52 A, a = 3.8 
A, and in-plane 87 and out-of-plane 88 conductivities yield 
£ofc ~ £c = 4.8 for the high- frequency dielectric con- 
stants, which determine the plasma frequency. For this 
choice in parameters we obtain £2 p ; = 0.914 eV at optimal 
doping in agreement with Ref. 87. In order to calculate 
the screened el-ph vertex the static dielectric constants 
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must be used in Eq. (24). Here, we take e c (0) = 10, 
which is obtained from optical measurements, 88 and set 
£ab(0) = 4e c (0). This choice in ratio results in a fac- 
tor two reduction in the breathing vertex at q = (0, 7r) 
which is consistent with Ref. 32. Additionally, we have 
set il — 60 meV in Fig. 6a. The real part of II^i rises 
as q 2 before abruptly falling as the imaginary part rises 
when the condition v F q = fi is satisfied. These kinematic 
constraints give similar forms for the mixed polarizabil- 
ities II fli i, although the detailed q dependence may be 
slightly different. 

In Fig. 6b-f several renormalized (screened) el-ph ver- 
tices g(k.AN,Qx,Qpl) are plotted for each of the modes 
considered in this work. Here k^^v = (t, 0) and the mo- 
mentum transfer is taken along the zone face q = {q x , 0) 
and g is plotted for several values of q z as indicated. 
To illustrate the effects of poor screening of the long- 
range Coulomb interaction, in Fig. 6b, a momentum- 
independent bare coupling g = 1 (arb. units) is also 
shown. For a momentum-independent vertex Tl g ,i — 
gHi,i. In this case, for 51 <§C £l p i(q) and q z = screen- 
ing is perfect, leaving only a small el-ph interaction at 
large q 2 £>- This can be seen as the black line, which cor- 
responds to q z = in Figs. 6b-f. For finite q z a small 
cone of wavevectors satisfies f2 > f2 p j(q) where screening 
is inoperable and the renormalized vertex recovers the 
bare value. This enhancement of the small q coupling 
occurs over a small range of q- vectors until (q 2 £>, Qz) sat- 
isfy flpi (q) = Qph , where a logarithmic divergence occurs 
which is cut-off by damping. For increasing q z this con- 
dition moves to a progressively larger q 2£) and a window 
of momentum points opens where the vertex can be well- 
represented by its bare value. 

Likewise, Fig. 6c-f plot the renormalized couplings 
for the Ai g , Bi g , apical and breathing modes, respec- 
tively. For large in-plane momentum transfers q 2D the 
c-axis couplings are largely unaffected by screening due 
to the fall-off of the polarizability at large <\ 2 d- How- 
ever, for the breathing branch, where the bare coupling 
weights large momentum transfers, the growth of the 
mixed polarizability leads to an overall suppression of 
the screened vertex at large momentum transfers. For 
small in-plane momentum transfers - relevant for c-axis 
Raman active phonons considered here - screening is in- 
effective and the coupling is anomalously enhanced over 
large q coupling. Again, for the breathing branch, only 
small effects are noticed at small q 2D due to the nature 
of the deformation-type coupling for these modes. Since 
we have neglected any fermionic dependence on the bare 
breathing couplings, the mixed polarizability is equal to 
the pure charge polarizability in this case. 

Using Eq. (22) the momentum-dependent el-ph cou- 
pling strength A„(q 2L >) 



K(<i3D) 



2N F 



k,q z 



(|g„(k,q,Sl p; )| 2 ) (25) 
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FIG. 7: (Color online) Plots of the momentum dependent 
coupling A„(q 2D ), Eq. (25), for v — A\ g (a), B\ g (b), apical 
(c) and bond-stretching (d) modes. The coupling are plotted 
on a logarithmic scale for q 2I) to highlight the region of poor 
screening. 



and bond-stretching modes shown in panels Fig. 7(a)- 
(d), respectively. One can clearly see that when summed 
over all out-of-plane momentum transfers q z , the net 
coupling is on the order of the bare coupling for small 
q 2 £i, while the coupling at large q is suppressed. The 
fermionic momentum dependence of the bare vertices 
for the Big case noticeably alters A(q 2£) ) for momen- 
tum transfers along the BZ diagonal where it is largely 
projected out. The A\ g and apical A(q 2D ) are quite sim- 
ilar even though the anisotropy of the bare interaction 
is substantially different and weights different regions of 
the BZ. For all the c-axis modes, coupling is small for 
Q2D = Qaf = {ir/a, 7r/o). Since we have neglected Hub- 
bard short-range Coulomb repulsion, which further sup- 
presses coupling at large momentum transfers, 35 Raman 
c-axis modes are thus not expected to appear in trans- 
port measurements. Likewise, the breathing modes are 
altered at small q 2D , but the overall coupling does not 
give much weight for these transfers. Since the apical 
and Big modes strongly favor coupling to the anti-nodal 
fermions, we may conclude that the underlying density 
of states, combined with the fermionic momentum de- 
pendence, strongly influences the magnitude of the cou- 
plings at small q 2D ; the phonon and eigenvectors control 
the actual angular dependence while the nature of the 
coupling - electrostatic or deformation - controls the |q|- 
dependence. 

This q-dependence strongly affects the d-wave pro- 
jected pair interaction as well. In Fig. 8 we plot the 
d-wave projected coupling 



is plotted in Fig. 7 as a function of q 2D (on a log scale to \ v d (q 2£) ) = 

highlight the small q behavior) for v = Ai g , Bi g , apical 
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FIG. 8: (Color online) Plots of the d-wave coupling \v,<i(<12d)i 
Eq. (26), for v — A\ g (a), B\ g (b), apical (c) and bond- 
stretching (d) modes. The coupling are plotted on a logarith- 
mic scale for c\ 2D to highlight the region of poor screening. 



for the same modes considered in Fig. 7. As with Fig. 
7, A ! y ! d(q 2r ,) has been plotted on a log scale in order to 
highlight behavior at small q 2D . One can see straightfor- 
wardly that the phonons that strongly favor small q scat- 
tering largely promote d-wave pairing (the c-axis modes) 
while phonons that favor large q-scattering, such as the 
bond-stretching modes, are detrimental. Furthermore, 
for the bond-stretching modes, while a large region of 
q space supports pairing, the large weight near q = Q 
dominates the coupling. The consequence of poor screen- 
ing for finite q z transfers, which accentuates small q cou- 
plings, is to enhance the overall d-wave coupling com- 
pared to q z — as well as standard Debye screening 
(flpi ^> J7 y ) for c-axis modes, and diminishes the repul- 
sive part for the bond-stretching modes. 

To visualize the effect of screening on the fermionic de- 
pendence of the coupling, we plot in Fig. 9 the screened 
X u (k FS ), derived from Eq. (19) where g has been re- 
placed with g from Eq. (22). Here, to mimic the effect 
of doping in the cuprates, we have varied the plasma 
frequency to smaller values in accordance with exper- 
iment, reflecting the increased insulating behavior and 
ionicity along the c-axis that occurs with underdoping in 
the phase diagram of hole- doped cuprates. The values 
for the total coupling A 2 = J2v(^v{k)) and A^ defined in 
Eq. (20) are shown in Fig. 10. 

Screening causes several noticeable effects compared to 
the unscreened case shown in Fig, 4. In order to mimic 
the effects of doping the screened el-ph vertex for the 
four branches is shown for £l p i ranging from 290 to 914 
meV. (The black arrows indicate the direction of increas- 
ing ilpi.) As the plasma frequency is increased the A lff , 
apical and breathing branches can be screened more effec- 
tively and the overall vertex is lowered around the Fermi 
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FIG. 9: (Color online) Plots of the screened e-ph A(1cf) for 
fermions at the Fermi level as a function of 9 for several values 
of the plasma frequency Q, p i. The value of £l p i varies from 
0.289 eV (black, solid) to 0.914 eV (purple, dashed). The 
black arrows indicate the direction of increasing tt p i. The 
parameters used to generate this hgure are defined in the text. 



surface. In the case of the A\ g and apical c-axis branches, 
the largest effect occurs in the anti-nodal region where 
the bare vertex is largest. In the case of the B\ g branch 
there is an anomalous anti-screening which occurs and 
the B\g vertex in the antinodal region is enhanced with 
increasing £l p i. This non- intuitive result is due to the 
out-of-phase oscillations of the oxygen modes since the 
only the difference between the bare vertices for the A ig 
and Big branches is the phase of the phonon eigenvectors. 
From these results it is evident that the self-energy due 
to coupling to c-axis phonons should redistribute weight 
around the Fermi surface as the number of doped holes 
varies. Thus, we infer that in the cuprates a window in q- 
space at small q 2 D occurs in which the el-ph interaction 
can still be quite large and avoid screening. This would 
not occur if the material was fully conducting along the 
c-axis and this effect increases with underdoping as fl p i 
is reduced with decreasing hole concentration. 

In contrast, since poor screening does not affect large 
momentum transfers, the strength of the breathing cou- 
pling dramatically increases for smaller J7 p j. This im- 
plies that the electron self-energy contribution from the 
breathing branches should grow with underdoping. Once 
again, however, we remark that this finding must be 
viewed with caution since large q-behavior is governed 
strongly by Coulomb interactions. 

Finally, in Fig. 10 we plot the total coupling A Zj re- 
sulting from the screened el-ph vertex as a function of 
fl p i. Following the trends in Fig. 9, the total coupling 
in both channels decreases with increasing fl p i for the 
in-plane breathing branch, as well as the apical and A\ g 
branches, while the coupling to the B^ g branch is en- 
hanced. In the case of the d-wave coupling, the attractive 
interaction of the A\ g branch is largely cancelled by the 
repulsive interaction of the breathing branch. Therefore, 
the total attractive interaction for d-wave pairing, which 
is primarily provided by the B\ g branch, is expected to 
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FIG. 10: (Color online) Plots of A z (left) and A^ (right) as a 
function of the plasma frequency Q p i. The parameters taken 
are denned in the text. 



be enhanced with progressive overdoping. 

The transition from insulating to more metallic behav- 
ior, reflected in Q p i, is one aspect of the doping depen- 
dence expected for el-ph coupling. However, this is not 
the only change expected to occur with doping. As a 
sample is doped away from half-filling the character of 
the pd-a* band is expected to change as spectral weight 
redistributes itself within the quasiparticle and lower and 
upper Hubbard bands. 89 This will be reflected in the pa- 
rameters Aq, A a and Aq u . Since these parameters enter 
the el-ph vertices in the fourth power, doping-induced 
changes in the bandstructure can have a large impact on 
the strength of the el-ph coupling. The value of the lo- 
cal crystal fields are also expected to vary with doping. 
This was shown in Ref. 90, which examined the role of 
an interstitial oxygen dopant in the SrO/BiO layers of 
Bi-2212 and found that the dopant locally enhanced E z 
at the planar oxygen site by up to a factor of five. The 
structural details of the crystal can also induce changes 
to the character of the band at Fermi level (Ao, Ac u 
and A a ), through variations in the Madelung energies of 
the atoms and the degree of hybridization between the 
CuC>2 plane and off-plane atoms such as the apical oxy- 
gen. Due to these considerations the overall strength of 
the el-ph coupling can have a complex dependence on 
both the carrier concentration of the CuC>2 plane as well 
as the composition and structure of the material. In the 
next section we examine the materials dependence of the 
el-ph interaction by conducting a systematic examination 
of the role of structure in determining the strength of the 
interaction across various families of high-T c cuprates. 



V. MADELUNG ENERGIES AND LOCAL 
FIELDS 

The environment around the CuC>2 plane influences 
electron dynamics in the plane via the Madelung ener- 



gies and local crystal fields. As a result, local symme- 
try breaking plays an important role in lattice dynam- 
ics. For example, X-ray measurements on YBa2Cu30y 91 
have shown that the CuC>2 planes are statically buckled, 
and a linear el-ph coupling results from the breaking of 
local mirror plane symmetry. LDA calculations for the 
infinite layer material CaCu02 also found evidence for 
static buckling of the plane due to the steric interactions 
among the oxygen 2p z orbitals which produces a sub- 
stantial linear B\ g el-ph coupling. 24 Among the different 
contributions to local symmetry breaking, this section 
is devoted to an investigation of the magnitude of the 
Madelung energies and local fields across families of the 
cuprates as a mechanism of local symmetry breaking. 

We remark that LDA investigations of el-ph coupling 
are limited mostly to stoichiometric YBCO, with virtual 
crystal extensions to LSCO and CaCuO. LDA bandstruc- 
ture, being metallic in nature, will tend to screen charge 
variations whereas those charge variations, particularly 
off the C11O2 plane, may be unscreened and result in 
strong local modifications of band parameters triggering 
nano-scale inhomogeneity as well as strong variations of 
el-ph coupling. This was demonstrated in Ref. 90, where 
it was shown that the presence of an interstitial oxy- 
gen dopant's unscreened charge in the SrO/BiO layers 
produced large local enhancements in the crystal field 
strength at the planar oxygen sites. This resulted in 
a local increase in the strength of the coupling to the 
A\ g and B\ g branches, which in turn reduced the energy 
gap and broadened the spectral features in the hole addi- 
tion/removal spectrum as well as locally increased J via 
the gain in lattice energy. 

Here, we focus our calculations in the ionic limit, where 
the ions are represented as point charges, and evaluate 
the electrostatic sums using Ewald's method. We show 
that the Madelung potential landscape is sensitive to 
the details of the unit cell, producing material-dependent 
variations in the field strength at the planar and apical 
oxygen sites. For the field at the oxygen site of the outer- 
most CuC>2 plane, which determines the coupling to the 
A\ g and B\ g modes, this variation mirrors the observed 
variations in T c suggesting a link between these quanti- 
ties. The use of formal valences in the electrostatic calcu- 
lations reflects the stociometric compounds and therefore 
the estimates obtained from their use corresponds to the 
parent compounds. However, in light of the findings of 
Ref. 90, we also address doping-induced changes in the 
crystal field strength in Bi-2212 in order to assess how 
these fields are expected to vary with doping. Finally, 
the variations in X z ^ that arise as a function of material 
due to the electrostatic considerations presented here are 
also discussed. Here we focus on the B lg mode, which 
has the largest contribution to d-wave pairing and take 
in into account variations in the local field, relative shifts 
in site energies, and changes in the atomic character of 
the band, all of which are directly affected by the details 
of the crystal structure. 
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FIG. 11: (Color online) The HgBa2Ca n _iCu„04n+a (n = 
1 — 5) unit cell showing the locations of the mirror plane (yel- 
low, MP), low symmetry breaking plane (green, LSB) and 
high symmetry breaking planes (red, HSB). The degree of 
symmetry breaking establishes the strength of the electric 
field which couples to the c-axis phonons. 



A. Formal Valances: The Parent Compounds 

A structural view of the lattice across the Hg family, 
as shown in Fig. 11. Starting with n = 1 shows that 
the coupling to planar oxygen c-axis phonons must be 
weak, while couplings involving the apical oxygen phonon 
can be strong due to the asymmetry along the c-axis 
at the apical site (Cu0 2 on one side and HgO on the 
other). For n > 1 the outer layer does not lie in a mirror- 
plane and the coupling strength depends on the local 
asymmetry (in this case BaO on one side, Ca on the 
other). However, the innermost plane may lie in a mirror 
plane when n is odd, while the intermediate planes will 
generally have a weaker degree of symmetry breaking in 
comparison to the outermost plane. Thus, in the context 
of el-ph coupling, the most strongly coupled modes will 
be those involving displacements of the apical and planar 
oxygen atoms in the outermost layers. In Hg compounds 
the local asymmetry around the outermost plane is due 
to Ca and BaO structures on opposite sides of the CuC"2 
plane. In YBCO it is due to Ca and SrO structures. 
Such changes in composition will affect the magnitude of 
the crystal field and these differences can be quantified 
by examining the Madelung potential for the apical and 
planar oxygen sites for the Hg, Tl and Bi cuprate families. 

The electric field strength can be regarded as the first 
derivative of the Madelung potential, which consists of a 
contribution of all electrostatic interactions between ions 
in the solid. For a test charge q t , located at position r, 
the Madelung potential $ is defined as: 



$ 



(Ji 



dt 



, pV) 



(27) 



where p is the charge density. In the ionic limit the charge 
density is comprised of point charges p(r) = J^ 8{r— Rj), 
where the sum i is taken over the atoms of the crystal 
located at R^. 

There are several common methods used to calculate 
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FIG. 12: (Color online) The Madelung potential at the Cu0 2 
planar oxygen site as function of distance from the plane of 
mirror symmetry: (a) HgBa2Ca n _iCu„02n+2 (n = 1-6, red 
circles), (b) Ti2Ba2Ca n _iCu n 02n+2, (n = 1-4, blue squares), 
(c) Bi2Sr2Ca n _iCu n 02n+2 (n — 1-3, green triangles), (d) An 
overlay of the Hg, Tl and Bi plots rescaled by the lattice 
spacing, a. The black line is a fitted curve of the form a<3? = 
,4 + Bcos(27rz/c). 



the Madelung potential. A general and powerful method 
was given by Ewald, 92,93 which is used here to deter- 
mine the Madelung energies for the planar and apical 
oxygen sites in the outermost CuC>2 planes in the unit 
cell for Bi (n=l-3), Tl (n=l-4) and Hg (n=l-6) cuprates. 
The atomic positions are obtained from the known struc- 
tural data given in Refs. 94-107. In all cases for- 
mal valence charges for the ions have been assumed. It 
should be noted that for the materials considered herein, 
the surfaces are non-polar, as opposed to YBa2Cu307 
(YBCO). 93 Therefore the surface termination layer is not 
critical for our calculations. 

In the case of YBCO, assuming that the surface is ter- 
minated between the BaO and CuO sublayers, a local 
field between 0.8 - 2.1 eV/A was obtained, arising from 
the asymmetry condition of Y 3+ and Ba 2+ on either 
side of the Cu02 plane. 93 This polar surface contribu- 
tion overwhelms the periodic contribution from Ewald's 
method. Without the polar contribution, the coupling 
for example to the Bi g mode is underestimated when 
compared to that obtained from periodic LDA methods. 
In Rcf. 93 it was shown how deviations away from for- 
mal valences can give much different values for the local 
crystal fields. We have also found that the local elec- 
tric fields are more sensitive than the Madelung energy 
to small structural distortions and deviations away from 
formal valences. Therefore we conclude that while the 
overall field values should be viewed as being approxi- 
mate, we assume that the variations across the different 
materials can be qualitatively compared. 

The Madelung energies of the planar and apical oxy- 
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FIG. 13: (Color online) The Madelung energy of the apical 
oxygen site as a function of distance from the plane of mir- 
ror symmetry: (a) HgBa2Ca n _iCu n 02n+2 (n = 1-6, red cir- 
cles), (b) Tl2Ba2Ca n _iCu„02n+2, (n = 1-4, blue squares), 
(c) Bi2Sr2Ca n _iCu n 02n+2 (n — 1-3, green triangles), (d) 
An overlay of the Hg, Tl and Bi plots in addition of points 
obtained for LSCO and LBCO (black open symbols). The 
dashed and solid black lines are fitted curves of the form 
a<!> — A + B cos(2-7r z / c) . The dashed curve is obtained when 
all of the data is included while the solid curve is obtained 
when the La data points are excluded. 



gen sites are plotted in Figs. 12 and 13, respectively, and 
plotted as a function of the distance between outermost 
CuC>2 plane (or apical oxygen) and the mirror plane lo- 
cated at the center of the unit cell. The calculated apical 
Madelung energies are in agreement with those obtained 
in Ref. 15 for similar structural data, while the Madelung 
energies for the planar sites are, to the best of our knowl- 
edge, presented here for the first time. One can see that 
the Madelung energies for both the planar and apical 
sites increase when the outermost Cu02 plane lies fur- 
ther from the center of the unit cell. This occurs both 
when the number of planes per unit cell is increased as 
well as when the layers are spaced closer to the charge 
reservoir due to the larger spacer ions between the Cu02 
layers. Scaling the Madelung energies by the in-plane lat- 
tice constant a collapses the data for each cuprate family 
onto a single curve, as shown in subpanel (d) of Figs. 
12 and 13. By symmetry, the electric field falls to zero 
at a mirror plane implying that $(z/c) must have zero 
slope at these points. If z is measured relative to the 
mirror plane, such as those indicated in Fig. 11, then 
the simplest function that satisfies these boundary con- 
ditions is a<t>(z/c) = A + B cos(2ir z / c) , and serves as an 
overall guide for describing the behavior of the data. For 
the planar oxygen sites we obtain A = 66.3, B = —24.0 
while for the apical oxygen sites A = 81.2, B = —23.1 
(A = 76.0, B = —32.2 when the points for the La family 
are excluded), in units of eV-A. 
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FIG. 14: (Color online) (al)-(a3) The magnitude of the elec- 
tric field at the apical oxygen site of the Hg, Tl and Bi families 
of cuprates, respectively. The data is plotted as a function of 
the apex distance from the plane of mirror symmetry located 
at the center of the CUO2 planes. (bl)-(b3) The magnitude 
of the electric field at the planar oxygen site of the outermost 
CUO2 plane for the set of same compounds. (cl)-(c3) The 
strength of the electric field scaled by the unit cell lattice pa- 
rameters for the planar and apical oxygen fields, respectively. 
The dashed lines represents a fit a 2 E(z/c)/c = Aam(2-irz/c), 
where A = 7.1 and 83.9 eV for the planar and apical sites, 
respectively. 



Fig. 14 shows the local E fields at the apical and planar 
oxygen sites for the outermost Cu02 plane in the unit cell 
for a number of cuprate families having different numbers 
of layers. For the single layer cuprates the CuC>2 plane 
lies in a mirror plane located at the center of the unit cell 
resulting in a zero field at the planar oxygen sites. Thus 
there is no el-ph first order coupling to the A lg and B\ g 
branches in the single layer materials. This is no longer 
the case for n > 1. Examining the largest data set, that of 
the Hg family (n = 1—6), the E'-field at the planar oxygen 
site rises from n = 1 and peaks for n = 3 where the 
location of the outer plane of the cell is z/c ~ 0.2. As n is 
increased beyond this point the magnitude of the crystal 
field begins to decrease. The data for the Bi- and Tl- 
families are consistent with this trend - generally the local 
fields rise and follow the associated Madelung energy as 
the outermost plane moves further from the center of the 
cell, and falls once the outermost plane moves closer to 
the edges of the unit cell where the Madelung energy 
profile flattens. Since the apical oxygen does not lie in a 
mirror plane for any cuprate, a finite field is allowed even 
for n = l. Behaviour similar to the planar oxygen site is 
also seen with a peak field strength generally occurring 
for n = 3 and decreasing beyond this point. 

The scatter in values is largely due to uncertainties 
in structural data as well as the oversimplified use of 
formal valences for all materials. However, the results 
are similar among the different cuprates if the position 
of the outermost plane relative to the center of the cell is 
considered. Figs. 14cl and 14c2 show that the local field 
may be collapsed onto a sinusoidal curve as a function of 
c-axis distance when rescaled to account for the unit cell 
volume variations. 
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The presented E-fic\d values are sensitive to the struc- 
tural data, particularly for the planar oxygen sites where 
the degree of static Cu-0 bond buckling is relatively un- 
certain. Since a static buckling further breaks mirror 
plane symmetry across the CuC>2 plane, small variations 
in the degree of buckling can have a large impact on the 
field values obtained while only having a minor effect on 
the Madelung energies. For example, neglecting buck- 
ling on the order of 0.001c in the n = 2 layer materials 
results in .E-field variations on the order of 20-40 per- 
cent for the planar oxygen sites, while the variations in 
the corresponding Madelung energies vary only by a few 
percent. We have found that both the Madelung energies 
and local fields for the apical oxygen are not as sensitive 
to the degree of static buckling in the CuC>2 plane. 

We now comment on the .E-field values reported here 
in the context of previous works involving coupling to the 
c-axis modes. Due to the fact that Ba 2+ and Sr 2+ lie on 
either side of the CuC>2 plane in Bi-2212, previous expec- 
tations were that the planar field in this system would be 
smaller than that found in YBCO and work by one of the 
authors on the Raman derived B\ g phonon lincshape con- 
firmed these expectations. 108 Since the T c 's of these two 
materials are roughly equivalent, it was concluded that 
the el-ph coupling was irrelevant to the pairing mech- 
anism. However, a re-examination of the data set, in 
view of the greater sample inhomogeneity and impurity 
effects in Bi-2212 compared to YBCO, indicates empiri- 
cally that the field may not be as small as previously con- 
sidered. Using the theory of Rcf. 62, the Fano line-shape 
of the Big phonon in optimally doped Bi-2212 108 can be 
fit equivalcntly with a value of the coupling, parameter- 
ized by A, which is an order of magnitude different if 
one assumes a large, intrinsic broadening of the phonon 
line. Since the values of A extracted from fitting Raman 
lineshapes are sensitive to disorder for values of A in the 
range of those found in YBa2Ca307_,5, a direct compar- 
ison with other compounds should again be considered 
qualitatively, with a view towards materials of compara- 
ble quality. Our Ewald calculations confirm that the field 
at the planar oxygen site for Bi-2212 (E z = 1.16 eV/A) 
is similar to YBa2Cu30y_,5. We show in the next section 
that this value is likely underestimated once considera- 
tions for doping are taken into account. 



B. Doping Dependence of the Electric Field 

We now briefly comment on the doping dependence 
of the local crystal fields. In discussing the trends in 
E z across the cuprates we have assumed formal valences 
for each lattice site, which is appropriate for undoped 
parent compounds. However, as was shown in Ref. 90, 
the doping process can significantly affect the strength of 
these fields resulting in a substantial local increase in A's. 
From a symmetry point of view, the movement of charge 
from the Cu02 plane to the charge reservoir layers with 
increased doping will enhance local symmetry breaking 



across plane. To model this effect in Bi-2212 we repeated 
the calculation for E z , but this time with the charge of 
the Cu02 planes uniformly raised by 0.15e, correspond- 
ing to optimal doping. In order to maintain charge neu- 
trality the charge of the BiO and SrO layers, where the O 
dopants are known to sit, 110 were uniformly reduced by 
the appropriate amount. Using this crude model we find 
the E z field at the planar site is raised to 3.56 eV/A while 
the field at the apical site is reduced to 16.33 cV/A, in 
agreement with the general trends previously reported. 90 
In section III we were interested in determining the val- 
ues of \z t< p at optimal doping and therefore these are the 
values of the crystal field strength which were used in 
section IIIB to obtain the unscreened A z ^(k) and total 
A z ,0's. 

Throughout this paper we have discussed a number of 
aspects of el-ph coupling in the cuprates which are af- 
fected by doping the system. The local -E-field strength, 
the ionicity of the compound, characterized by Qpi, and 
the orbital character of the band at the Fermi level all 
modify the el-ph couplings and arc all affected by doping. 
Therefore, el-ph coupling in the cuprates is expected to 
exhibit a strong doping dependence. We also note that 
many of these parameters are also expected to vary be- 
tween materials comprising the various families and we 
turn our attention to these trends in the following sec- 
tion. 



C. Material Dependence of \ z ^ for the B\ g branch 

We end this section by examining the systematic vari- 
ation in the el-ph coupling strength expected from ma- 
terial to material. Since our goal is to understand the 
materials variations in T c we focus our attention on cou- 
pling to the Big branch, which provides the largest con- 
tribution to pairing and whose coupling is the strongest 
after the effects of screening are taken into account. As 
discussed previously, there is a empirical relationship be- 
tween T c and the distance of the apical oxygen from the 
Cu02 plane and this relationship this has generally been 
tied to the effective increase in the next nearest neighbour 
hopping t' through hybridization effects associated with 
the Cu 4s orbital. These hybridization effects will also 
alter the character of the pd-a* band crossing the Fermi 
level and can therefore affect the overall strength of the 
coupling to all of the phonon branches considered in this 
work. For the Bi g branch, such changes are reflected in 
Ao and, since this enters the total coupling as the fourth 
power, changes in hybridization can have a large affect 
on the total coupling strength. 

In Fig. 15a the systematic variation of Ao is presented 
as a function of T c for the same materials considered 
in the Fig. 14. To estimate Ao we consider the same 
five-band model used in section IIIB, but assign the site 
energies based on the Madelung calculations presented 
here and the energy level scheme of Ref. 15. Using hole 
language, we set e^ = 0, e p — A, e a = A + A$^/e(oo) 
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FIG. 15: (Color online) (a) The planar oxygen character at 
the Fermi surface Aq as a function of the material's T c . (b) 
The corresponding local _E-field at the planar oxygen site of 
the outermost Cu02 plane in the parent compound, (c) The 
value of \ z for the B\ g modes of the outermost plane for each 
material. 



and e s = e d — 7 . Here, A is the charge transfer energy, 
which is related to the difference in Madelung energies 
between the Cu and planar O sites A$j\/ = <I>o — 3>c« 
and is given by 



A = 



A$ 



M 



e(oo) 



Icu{2) + A {2) 



(28) 



where Ic«.(2) and Ao{2) are the second ionization and 
second electron affinity energies for the Cu and O sites, 
respectively. The factor of e 2 /d p represents the contribu- 
tion of the Coulomb interaction between the introduced 
electron-hole pair. Following Ref. 15, we take the dielec- 
tric constant e(oo) = 3.5 and Icu(2) — Ao(2) + e 2 jd p = 
10.9 eV. In determining e a , A§a = &A — ®o denotes the 
Madelung energy difference between the apical and pla- 
nar oxygen sites. In addition to varying the site energies 
between materials, the values of the hoppings t' and t 



ps 



are also adjusted using the prescription of Ref. 109. This 
allows us to account for changes in hybridization due to 
the different apical distance given in the experimental 



structural data. Furthermore, we neglect any variation 
of the a and b lattice constants and hold the in-plane hy- 
bridizations fixed throughout. Finally, correlations are 
handled at the mean field level. The overall effect of the 
mean field corrections is to shift the site energies 






U dd (nf)/2 + ±U pd (n?) + U pd (n?) 



U, 



Mi \ 



^)/2- 



2U, 



pd\ 



(29) 



= e a + U pp (n^/2 + U pd (nf) 



and e s is fixed to maintain the difference e d — e s — 7 eV. 
The Hubbard repulsions are set to canonical values with 
(in eV) U dd = 8, U pp = 4 and U pd = 1. In Eq. (29), 
nf denotes the total number operator for orbital a and 
site i and the paramagnetic solution for (nf) has been 
assumed. 

The results for A<j, E z and \ z ,4> are shown in Fig. 
15. The E-fields presented here are the same as those 
used in Fig. 14, and were obtained by assigning formal 
valences for each ion. These field values reflect best the 
crystal field of the undoped parent compounds and there- 
fore Fig. 15 does not include the effects of doping. As 
we have already seen, even in this simple model, the re- 
distribution of charge from the charge reservoir to the 
Cu02 plane can alter the crystal field values resulting in 
large changes in the value of A z>( j,. In the case of Bi-2212 
this produced an order of magnitude change in the value 
of A Zj for the factor 4 change in E z in the doped lat- 
tice. Furthermore, the ionic point charge model neglects 
any covalent nature of the bonds which will smooth the 
charge distributions. Therefore, the field strengths in real 
materials may be quite different from the values reported 
here, especially in the doped systems. However, by us- 
ing the same method for each material we can obtain 
valuable comparative information about the crystal field 
variations between materials. 

As shown in Fig. 15a, there is a correlation between 
the material's maximum T c at optimal doping (apical 
distance) and the planar oxygen character Ao- Again, 
this can be understood in terms of the increase in the 
effective planar 0-0 hopping via the Cu 4s as the api- 
cal oxygen is further removed from the Cu02 plane. A 
similar correlation, shown in Fig. 15b, occurs between 
E z at the planar oxygen site and T c for the reasons pre- 
viously discussed. Since the total couplings \ z ^ for the 
Big branch are proportional to the product A^E Z a clear 
correlation between the material's T c and the B\ g cou- 
pling naturally arises. This is shown in Fig. 15c. 

The correlation between T c and the strength of the 
coupling to the B\ g branch, which provides the largest 
contribution to d-wave pairing, provides a natural frame- 
work for thinking about materials variations in T c . Al- 
though it is clear that the overall strength of the el-ph 
interaction is too small to support HTSC on its own, 64-66 
phonons are not precluded from contributing to pairing, 
enhancing the pairing correlations from another domi- 
nant interaction. Such a possibility has been considered 
in previous works that have examined the contribution 
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of el-ph coupling in conjunction with antiferromagnetic 
spin fluctuations. 111,112 In such a scenario, a baseline T c 
is set by the dominant pairing mechanism, typically iden- 
tified with antiferromagnetic spin fluctuations. Such an 
interaction is likely to be governed by the properties of 
the CuC>2 plane, such as the charge transfer energy A 
or antiferromagnetic exchange energy J, and therefore 
the strength of this interaction is likely independent of 
the material. T c is then further enhanced by the weaker 
contribution from the el-ph interaction which, as we have 
shown, is strongly materials dependent. This picture nat- 
urally explains the variations in T c that are observed 
across materials and provides a microscopic connection 
between the structure and composition of the unit cell 
and the T c of the material. 

In the next section wc continue examine such a multi- 
channel model and present simplified considerations for 
T c in order to assess the degree to which phonons can 
enhance a dominant interaction at high energy. As we 
will show, moderate values of A Zj 0, such as those given in 
Fig. 15c, can lead to sizeable enhancement to T c which 
is greater than the T c that would result from phonons 
alone. Phonons are therefore capable of providing a suf- 
ficiently large enhancement to T c despite the relatively 
low value of their total coupling \ z ,<p. 



VI. CONSIDERATIONS OF MULTIPLE 

PAIRING CHANNELS: T c AND THE ISOTOPE 

EFFECT 

In this section we turn our attention to considerations 
for T c by considering contributions to pairing from mul- 



tiple channels. T c , one of the most difficult quantity to 
calculate correctly, varies across the cuprate families by 
two orders of magnitude even though the same Cu02 
building block is present. However, as we have demon- 
strated, the total el-ph coupling for the B\ g branch shows 
a strong materials dependence which is empirically cor- 
related with the observed value of T c at optimal doping. 
Therefore it is of interest to explore how el-ph interac- 
tions may provide a material dependence to T c arising 
from the material conditions along the c-axis and poor 
screening. 



To calculate T c we consider a multi-channel boson 
exchange model for pairing, which is a straightforward 
generalization of the Eliashberg equations. A similar 
approach has recently been used 113 to account for the 
qualitative structures in the phonon-modulated density 
of states of the cuprates as observed by STM. 45 Implicit 
in this calculation is the assumption that the dominant 
pairing channel can be described in terms of a boson 
exchange mechanism and such a scenario is the sim- 
plest one which can be adopted in examining the idea of 
phonon-enhanccd HTSC. However, such a model can be 
perhaps be justified based on several recent works. For 
example, studies of the penetration depth in electron- 
doped cuprates using a model with coexisting AFM 
and SC order has produced good agreement with ex- 
periment indicating the importance of AFM in these 
systems. 114 Furthermore, a DMFT study of optical con- 
ductivity in a moderately correlated single-band Hub- 
bard model, 115 has concluded that AFM with moderate 
correlations is the correct picture for describing the un- 
derdoped cuprates rather than strongly correlated Mott 
physics. Finally, a study of the t-J and Hubbard mod- 
els has shown that the d-wave pairing in these models 
have a significant contribution from a retard component 
with an energy set by a few times J or a fraction of 
U in the t-J and Hubbard models, respectively 116 This 
result indicates that the pairing interaction in these mod- 
els is largely retarded in nature with an energy scale set 
by the dynamics of the magnetic excitation spectrum. 
These theoretical studies lend some support to the spin- 
fluctuation mediated picture adopted here although the 
details may differ. Nevertheless, this model provides a 
simple means to evaluate the role of phonons in modify- 
ing T c in a multi-channel model. 



For this calculation we assume that each channel can 
be described in terms of a boson exchange mechanism. 
Individual channels are indexed by v and the el-boson 
contributions to the mass renormalization Z(uj) and 
anomalous self-energy </>(#, w) are parameterized by X VyZ 
and A„,0, respectively. Making the standard approxima- 
tions of a structureless band and cylindrical Fermi sur- 
face, the multi-channel Eliashberg equations are 
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A(oj)Z(u) = ]T 



\ v ,4,n v , dJ 



v /cj' 2 -A 2 (w)cos 2 (20) 
A(w')cos 2 (26») 



yV 2 - A 2 (a/) cos 2 (2(9) 




Here we have assumed a separable form for the anomalous self-energy 4>(9,oj) = (j>(u>) cos(2#) and A(u>) 
is the complex gap function. The kernels in Eq. (30) are given by 



(30) 

<P(lo)/Z(oj) 



K^(w',u) = [f(-w')+n(n v )] 



1 



uj + f2„ + u>' 



± 



cj — fi^ — u 1 



+ [f(uj')+n(n„)} 



■-SK 



± 



fl v — i 



where / and n are the Fermi and Bose occupation factors, respectively. Following Rcf. Ill we adopt a double 
square-well model for the frequency dependence of A(uj) and Z(ui), as shown in Fig. 16. Here, we identify v — 1 
with the f2i ~ 35 meV B\ g phonon branch, and v = 2 is identified with antiferromagnetic spin fluctuations with 
fi 2 = 2 J ~ 260 meV, reflecting the top of the magnon band in the case of a long-range antifcrromagent. Next, wc 
follow a standard set of approximations 117 and obtain the coupled equations 



/A 
dJ — * tanh 



for fl\ > w > and 



2fc fc T c 



/■ A 2 

+ (Ai,0 + A 2 ^) / dw — tanh 
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ctcj — tanh 



(32) 



for f2 2 > w > f2i. Here Z\ = 1 + Ai )2 + A 2 . z and Z 2 = 1 + A 2 . 



r 



Eqs. (32) and (31) are then solved for the non-trivial 
solution for Aj and the resulting expression for T c is 



k B T c = 1.134M^~ 



7 f2 2 exp 



1 + A M + A 2|Z 

Al,<A + A 2 ,<A 



where 
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1 + Ai, 



1 + A, 



Ai,0 + A 



2.ffl 



(33) 



(34) 



This expression recovers the McMillan result (we have 
neglected /i*) for the single well case in cither limit 
when Xi. z ,Xi^ or A 2iZ ,A 2i are set to zero. Other ex- 
pressions have been obtained for T c using different forms 
of the two-well model and neglecting the difference in 
contributions from each of the momentum channels (i.e. 
A z = A^) while taking the spin-fluctuation contribution 
to be repulsive. 118 ~ 120 We refer the reader to Ref. 120 
for a more thorough discussion of the various proposed 
expressions. 

In order to assess the phonon contribution to T c , we 
now assume that spin-fluctuations provide the dominant 
source for d-wavc pairing and parameterize this mode 



T c = 62 K in the absence of phonons. For the phonon 
mode we assume fix = 36 meV and hold Ai. z = 2Ai ) ^. 
In Fig. 17a the resulting T c as a function of el-ph cou- 
pling is presented. For reference, the value of T c obtained 
from phonons alone is also shown, multiplied by a factor 
of 10. It is important to note that even small el-ph cou- 
pling can enhance T c considerably over the single-channel 
spin fluctuation model. Furthermore, the degree of en- 
hancement depends on the ratio of A 2i< ^,/A 2)2 assumed for 
the spin fluctuations. In Fig. 17b we plot T c for the 
two-mode model while varying ratios of A 2j 0/A 2iZ . In 
each case we have adjusted A 2i such that the baseline 
of T c — 62 K is maintained for Ai^ = 0. In general, as 
A 2 .^/A 2iZ increases in the dominant channel, the overall 
gain in T c mediated by phonons also increases but with 
diminishing returns for larger A 2 .^/A 2 , z . An isotropic re- 
pulsion in the spin fluctuation channel is also expected 
to produce a similar effect. We also note that our results 
are limited by the questionable applicability of Migdal- 
Eliashberg theory in the cuprates and therefore T c itself 
is a poor metric of a theory to be compared to exper- 
iment. Nevertheless, as this enhancement is generic to 
multi-channel couplings, it may serve as a guide for un- 
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FIG. 17: The transition temperature T c as a function of 
el-ph coupling in the d-wave channel Ai,0, calculated in the 
two-square well model, defined in the text. The dashed line 
is T c expected for the phonon alone. 



derstanding the material dependence of T c itself. 

In Fig. 18 the isotope exponent a = —dlnT c /d\uM, 
where M is the mass of the oxygen ion, is plotted using 
the same parameters as in Fig. 17 with a = (1— 7)/2, and 
7 defined by Eq. (34), reducing the overall magnitude of 
a from the value for phonons alone. A similar reduction 
in a and enhancement of T c has been reported in Ref. 
112 when the phonon contribution to pairing was taken 
to be attractive in the d x i_ y i channel. Fig. 18 shows 
that the overall isotope exponent is small even though 
the phonon mediated enhancement to T c is sizeable and 
demonstrates that small isotope exponents are possible 
despite large enhancements in T c due to the el-ph inter- 
action. 

Some degree of caution is warranted in making a di- 
rect comparison to experiment where the isotope ex- 
ponent is small at optimal doping and increases with 
underdoping, 21 while our results show that the overall 
contribution of el-ph coupling to <i-wave pairing is ex- 
pected to increase with doping due to screening effects. 
Therefore one might expect that the isotope exponent 
should increase with doping, contrary to what has been 
observed experimentally. However, doping dependent 
changes are also likely to occur for a spin-fluctuation me- 
diated pairing mechanism. Since it is the relative contri- 
butions from each channel that sets the value of a a direct 
calculation of a expected experimentally is not possible 
until the details of the dominate pairing interaction are 
understood. Nevertheless, these calculations show that 
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phonons with relatively low el-ph coupling strengths can, 
in principle, play a substantial role in determining T c 
while producing small signatures in traditional metrics 
such as the isotope exponent. 

We end with a few comments about these results in re- 
lation to prior works. As we have already noted, the spe- 
cific results presented are limited by the applicability of 
Eliashberg theory and the Fermi liquid description of the 
cuprates. Such a picture may not be completely applica- 
ble, especially in the underdoped region of the phase dia- 
gram where the role of correlations becomes more promi- 
nent. As the role of correlations and its interplay with 
the el-ph interaction is better understood the details of 
the picture presented here may change. However, the 
phonon-mediated enhancement of T c is expected remain 
as a general phenomenon. 

We also note that the inclusion of correlations may 
further enhance the phonon contribution to T c through 
possible enhancements of the el-ph vertex or proxim- 
ity to AFM. For example, Recent DMFT work examin- 
ing phonon-mediated pairing in the fullerides 121 ' 122 has 
found that a phonon-mediated pairing mechanism, which 
produces an exponentially small T c when considered on 
its own, will be significantly enhanced in proximity to 
a metal-insulator transition. The mechanism via which 
this occurs is a density of states enhancement due to a 
narrowing of the quasi-particle band and Mott physics. 
This is a different origin than the enhancement reported 
here where we are concerned with the overall effects of 
coherence and band width of the paired electrons. Some 
aspects of the Mott physics is included in our considera- 
tions in the form of the breakdown of metallic screening 
along the c-axis (a Coulomb effect) however, we have not 
included the breakdown of carrier properties in the Cu02 
plane. This is an intriguing possibility but it is beyond 
the scope of our paper and we leave it for future work. 

Finally, it has recently shown that high values of T c 
are also obtainable from el-ph coupling without invoking 
a second dominant pairing interaction by departing from 
the conventional Fermi liquid framework. In Ref. 123 it 
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was shown that T c 's on the order of 100 K are attainable 
for A Z] values similar to those obtained here, provided 
the normal state behaves as a non-Fermi-liquid quantum 
critical metal. Although the simple approach outlined 
here does not hold in this scenario, much of the doping 
and materials-dependent variations in coupling still do. 



VII. SUMMARY AND CONCLUSIONS 

Although there appears to be a great deal of universal- 
ity in the hole-doped cuprates, large changes in T c can 
be achieved through changes in the chemical composi- 
tion and structure of the unit cell. As we have shown 
here, these changes are reflected in a number of proper- 
ties including the orbital character of the band crossing 
the Fermi level, the strength of local electric fields arising 
from structural-induced symmetry breaking, doping de- 
pendent changes in the underlying band structure, and 
ionicity of the crystal governing its ability to screen c- 
axis perturbations. Given the sensitivity of T c to the 
structural details of the crystal it is clear that the un- 
derlying mechanism(s) for HTSC must incorporate these 
elements in some way. In this work we have examined 
a number of aspects of coupling to oxygen modes and 
demonstrated that the overall coupling to these modes 
is quite sensitive to these factors. Therefore, the inclu- 
sion of el-ph coupling to these modes provides a natural 
means of linking the electronic properties of the Cu02 
plane to the structural elements of the material lying off 
plane and the carrier concentration. This picture calls 
for oxygen phonons to play a role in HTSC. 

In extending the previous works by some of the authors 
on el-ph coupling in the cuprates, we have formulated a 
theory for poor screening in these materials. Due to the 
quasi-2D nature of the cuprates with poor conductivity 
along the c-axis, we have shown that a window of small 
q opens in which screening is inoperable. This results 
in an overall enhancement of the total (i-wave projected 
coupling which enhances the phonon's ability to medi- 
ate d-wave pairing. With progressive doping, the total 
phonon contribution to A^ is enhanced while the total 
contribution from A z , which suppresses d-wave pairing, 
is screened away. As a result, doping the system away 
from half-filling results in an increased contribution to 
d-wave superconductivity mediated by el-ph coupling. 

A systematic examination of the variation in el-ph cou- 
pling across the cuprate families was also performed. Us- 
ing experimental structural data we determined the ma- 
terials variation in the total electric field strength and or- 
bital character of the pd-a* band at the Fermi level and 
the expected variations in el-ph coupling strength to the 
Big branch. A direct correlation between the strength of 
the coupling and the optimal T c for each material was 
observed. As a result, the materials variation in T c can 
be naturally accounted for through the inclusion of el-ph 
coupling to the presently accepted intrinsic planar mod- 
els. 



To be clear, we do not suggest that el-ph coupling 
alone accounts for HTSC as both the bare and screened 
values for the el-ph coupling strength arc insufficient to 
produce large T c 's on their own. Instead, we propose 
that the phonons work in conjunction with the presently 
unidentified dominant mechanism for HTSC. The sim- 
plified multi-channel model for T c presented here has 
demonstrated that moderate el-ph coupling, when com- 
bined with spin fluctuations, can produce enhancements 
to the total T c in excess of the contribution produced 
by phonons alone. Furthermore, this enhancement is ex- 
pected to be a generic phenomena and independent of the 
presently second mechanism. Therefore, the inclusion of 
the el-ph interaction can account for the large variations 
in T c observed across the cuprate as a function of mate- 
rial and doping. 

There are some open issues regarding el-ph coupling 
in the cuprates. First of all, our examination of the field 
strength as a function of material neglected any doping- 
induced changes to the field strengths. Our simplified 
treatment of the bi-layer system Bi-2212 showed that 
these effects can be substantial. Therefore, this problem 
requires a systematic examination of the doping process 
across the families of cuprates, perhaps with more so- 
phisticated approaches. This issue is complicated in the 
multi-layer systems with n > 2 where there is the pos- 
sibility for inequivalent dopings for the Cu02 planes. 49 
This can further modify charge distributions altering the 
strength of the crystal fields in each plane. 55 Any system- 
atic treatment of the doping induced fields should take 
this into account. There is also the additional question 
of the role of covalency in determining the overall field 
strength. 

Finally, we have investigated the role of el-ph cou- 
pling in the cuprates, taking into account the long range 
Coulomb interaction and its role in screening. However, 
the role of the short-range Coulomb interaction and its af- 
fect on the el-ph interaction remains an open and intrigu- 
ing question. We believe that the findings presented in 
this paper serve to further highlight the need for progress 
in this area. 
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Appendix A: Modifications of the Zhang-Rice 
Singlet 

Throughout this work the role of electronic correla- 
tions was neglected in deriving the anisotropic el-ph cou- 
plings arising from different charge transfer processes. In 
this appendix an analysis via a different route is under- 
taken by examining directly the role played by corre- 
lations. Previous studies on Hubbard clusters have il- 
lustrated the usefulness of exact diagonalization for ex- 
ploring low energy excitations, such as the Zhang-Rice 
singlet (ZRS), 124 as a function of hole number. In par- 
ticular, these studies have demonstrated a strong depen- 
dence of the ZRS energy and mobility, as well as the 
magnetic exchange J, on the planar charge transfer en- 
ergy A = e p — £d- 125 Ref. 15 pointed out that the sta- 
bility of the ZRS increases if the apical atom lies further 
away from the CuC>2 plane, which reduces the mixing 
of the ZRS with other low energy states of the cluster. 
Since the phonons under consideration directly modulate 
A and the apical Madelung energy, strong modifications 
should be seen in the ZRS and magnetic exchange. The 
cell perturbation approach in Ref. 126, which included 
in-plane oxygen breathing and buckling modes, showed 
that an effective simple Hubbard-like Hamiltonian can be 
deduced with renormalized parameters. 

Here, exact diagonalization studies of multi-band Hub- 
bard clusters are used to investigate how the low energy 
sector is modified by local, static atomic displacements 
in a pattern given by various phonon eigenvectors. The 
modifications of the ZRS and J due to el-ph coupling 
receive close attention in order to derive renormalized 
parameters for single-band Hubbard to t-J model Hamil- 
tonians. 

The clusters show in Fig. 19 contain two sets of or- 
bitals at the Cu site: Zd x 2_ y 2 and 4s. By symmetry, the 



hopping amplitude from the apical oxygen site to the Cu 
3d x 2_ y 2, t a d vanishes, while t as is non-zero. Thus, a gen- 
eral Hubbard model can be written as H = H^ in + H pot 
where the kinetic piece is 



Hkin = t pd y^ Qs[d\ a Pi,s,a + h.c] 



tap^2 Qs[a\.aPi,S,a + h. 



i,5,cr 



t as '^2 l [a\,<7 S i,<y + h.c] 

tps ^2 Qs[4,aPi,S.a + h.C.' 



(Al) 



with phase factors Q± x — —Q± y = ^1, Q'± x ,±y = 1 = 
Q'± XTy , Q± x = Q'j. y = ±1, and the potential piece is 
given by 



i.S 



H po t = UddJ^nitKi + UppJ^^l 

i 

+ td^fifv + e^hl 

+ ^(e + tf -£aKa 

+ ^2{^v + Ui,5-E p )nl 



S,^ n i,S,^ 



(A2) 



Here, a, ps, d, s (a\p\, d\ s') annihilate (create) a hole 
on apical oxygen, in-plane oxygen (S = x,y), copper d- 
and copper s-orbital, respectively, U is the displacement 
of a given ion and E is the associated local electronic 
field. This type of coupling is different than that consid- 
ered in previous studies, where the buckling modes were 
represented as an electrostatic modulation of the Cu site 
energies. 126 Since the 4s orbitals are rather extended we 
neglect the on-site Coulomb repulsion as well as inter- 
orbital interactions with the 3d x 2_ y 2 orbitals. 

There arc a total of 13 and 14 orbitals for the Cu20g 
and CU2O10 clusters, respectively. In both cases, the 
clusters contain seven holes in order to investigate the 
formation and derealization of a ZRS to a nearest neigh- 
bour copper-oxide plaquette (4 holes being almost ex- 
clusively located in the Cu 4s orbitals). Exact diago- 
nalizations in the S* ot = 1/2 sector, which contains the 
ground state, are performed since the Hamiltonian con- 
serves Sf ot . Parameters similar to those of Ref. 15 for 
the Hubbard model are used, with values representative 
of the electric field for the planar and apical oxygen sites 
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FIG. 20: (Color online) The effect of different modes on 
various Zhang- Rice parameters for the CU2O9 cluster: plots 
(a) and (b) show the effect of a given mode on the ground- 
state energy and on the effective ZRS hopping t with At = 
\t(5) - t(8 = 0)| and AE = E(6) - E(S = 0). (a) The effect 
of the A lg modes, (b) The effect of the B lg mode, (c) The 
effect of the apical mode on 2t. 



in bilayer cuprates (in eV and eV/Afor the fields): 



tpd — 1-13 


t — 2 

L ps ~~~ ^ 


DT) VJ . t: J? 


tap — u.zy 


/ — 2 


e s = -7 


€p € a z.y 


Upp = 4.1 


U dd = 8.5 


E p = 1.6 



E a = 16 

The diagonalization of the CU2O9 cluster yields the ef- 
fective hopping parameter t of the ZRS as well as the 
magnetic exchange J: the splitting between the ground- 
state and first excited state gives 2t, whereas the singlet- 
triplet splitting of the 6-particle problem gives J. The 
parameters t' and J' are analogously defined from the 
Cu 2 Oio cluster. In the absence of el-ph coupling we ob- 
tain the following values: |2t| = 0.63 eV, \2t'\ = 0.28 eV, 
J = 0.17 eV, and J' = 15 meV. 

Fig. 20 displays the effect of different phonon modes 
on the ground state energy, t and J as a function of 
static displacement 5. Although only static atomic dis- 
placements have been investigated in configurations cor- 
responding to the phonon eigenvectors, the differences 
between the relevant time scales associated with the elec- 
tron hopping and atomic motion justify such an analysis 
for the purpose of viewing a snap-shot of the modifica- 
tions of the ZRS parameters. (We also note that Ref. 
90 considered a similar calculation with fully quantized 
atomic motion for the B\ g and A\ g phonons at q = 
and reported similar results for J.) One clearly sees, 
for small displacements, the main effect on t is due to 
the Aig and B\ g modes. For the B\ g mode, the effective 




FIG. 21: (Color online) Variation of the effective Zhang- 
Rice hopping t' and magnetic exchange J' to next-nearest 
neighbor. The calculation performed for the CU2O10 cluster 
for Ai g -, Big- and apex-modes. 



hopping of the ZRS is modulated to first order in the dis- 
placement, while the energy is only modified at second 
order, in agreement with perturbative analysis. For the 
Ai g mode the first-order correction to the energy does 
not vanish and both AE and At oc S. The magnitude 
of the modulation in t is comparable to the B\ g mode 
with the same local field. The effect of the apical mode 
on t is very weak for S > 0. This is not surprising since 
the ZRS has little apical character if the orbital lies away 
from the CuC>2 plane. For 6 < the hybridization of the 
apical-p and Cu-d3 Z 2_ r 2 orbitals destabilizes the ZRS and 
thus the effective hopping t drops as in an earlier work. 15 
Also, both the ground state energy E and the effective 
exchange J are not affected much by this mode. 

It is clear, based on the calculations performed, that 
the three modes modify the effective hopping of the ZRS. 
However, the mechanism by which the apical mode af- 
fects the hopping is very different since the formation of 
the ZRS no longer takes place in the Cu02 plane if the 
displacement of the apical oxygen is towards the plane. 
The effective next-nearest neighbour quantities t' and J' 
(see Fig. 21) are affected less by the two planar modes. 

Summarized in Figs. 20 and 21, our results indicate 
that for both the A\ g and B\ g phonons, At as well as A J 
corrections have a linear dependence on displacement of 
about the same magnitude. For the apical phonon the 
magnitude is different due to the destabilization of the 
ZRS. The on-site energy is modified mainly by the A\ g 
ep-ph coupling; which can be quantified by looking at 
the corresponding modification of the ground-state en- 
ergy. Modifications to parameters AJ' and At' indicate 
again that the dominant contribution comes from the A\ g 
mode. However, this conclusion is based on a naive clus- 
ter approach and should be corroborated by other means 
of investigation. 

In terms of treatments of effective single-band models, 



2(> 



such as Hubbard or t- J, incorporating phonons, we note 
that our results indicate that the effective couplings for 
these modes are neither purely Holstcin-like nor can they 
be treated as bond phonons simply modifying the effec- 
tive ZRS hopping t: the effective parameters, such as t 
and J, and as a consequence, the Mott-Hubbard effective 
interaction U in a single-band approach, are all modi- 
fied. The actual modification of the effective U for static 



displacements corresponding to these phonon modes re- 
quires an investigation of different three-band clusters 
containing Cu d s , d 9 and d 10 and oxygen ligand con- 
figurations. This points out the complications involved 
in ascribing real phonon modes to any corresponding 
model phonons in a single-band approach, and illustrates 
the important interplay between el-ph coupling and elec- 
tronic correlations. 
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